(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     73453,       1527]
NotebookOptionsPosition[     72757,       1502]
NotebookOutlinePosition[     73115,       1518]
CellTagsIndexPosition[     73072,       1515]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"\n", 
  RowBox[{"(*", 
   RowBox[{"reading", " ", "in", " ", "the", " ", "human", " ", "data"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"LineNumber", "=", "1"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"IndexAllPre", "=", 
     RowBox[{"Import", "[", 
      RowBox[{
      "\"\</home/shu/workspace/Research/AMBER_Output_Test/data/IndexAll.mat\>\
\"", ",", "\"\<MAT\>\""}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"IndexAll", "=", 
     RowBox[{"Join", "[", 
      RowBox[{"First", "[", "IndexAllPre", "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"HipPosIndex", "=", 
     RowBox[{"Round", "[", 
      RowBox[{"IndexAll", "[", 
       RowBox[{"[", 
        RowBox[{"LineNumber", ",", "1"}], "]"}], "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"NSslopeIndex", "=", 
     RowBox[{"Round", "[", 
      RowBox[{"IndexAll", "[", 
       RowBox[{"[", 
        RowBox[{"LineNumber", ",", "3"}], "]"}], "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"TorsoIndex", "=", 
     RowBox[{"Round", "[", 
      RowBox[{"IndexAll", "[", 
       RowBox[{"[", 
        RowBox[{"LineNumber", ",", "3"}], "]"}], "]"}], "]"}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5206980829927273`*^9, 3.520698089238552*^9}, {
   3.523978947801125*^9, 3.5239789588067827`*^9}, {3.52397911342336*^9, 
   3.523979114631843*^9}, {3.5243169653528748`*^9, 3.524316973652472*^9}, {
   3.5243319180886717`*^9, 3.52433193045708*^9}, 3.524332052043282*^9, {
   3.5243336818168917`*^9, 3.524333682086383*^9}, {3.527031289451779*^9, 
   3.52703129229949*^9}, {3.527080033657362*^9, 3.52708006144129*^9}, {
   3.527080133566532*^9, 3.52708014024951*^9}, {3.5271009843940268`*^9, 
   3.5271010408501663`*^9}, 3.527250709910329*^9, {3.5272507404573*^9, 
   3.5272507718199787`*^9}, 3.52725082256019*^9, {3.527279016037776*^9, 
   3.5272790190211163`*^9}, 3.5272794890519753`*^9, {3.5272832222939243`*^9, 
   3.527283234403532*^9}, {3.5274518483084097`*^9, 3.527451848605681*^9}, {
   3.529074757490752*^9, 3.52907475844724*^9}, {3.529855356841095*^9, 
   3.529855356934772*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "choose", " ", "for", " ", "mean", " ", "data", " ", "or", " ", "fred", 
    " ", "data"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"DataType", " ", "=", " ", "\"\<mean\>\""}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "choose", " ", "the", " ", "type", " ", "of", " ", "hip", " ", 
      "position", " ", 
      RowBox[{"output", ":", "\[IndentingNewLine]", 
       RowBox[{"1.", " ", "hip", " ", "position"}]}]}], ";", "\n", 
     RowBox[{"2.", " ", "linearized", " ", "hip", " ", "position"}], ";", 
     "\n", 
     RowBox[{
     "3.", " ", "Dimensionless", " ", "linearized", " ", "hip", " ", 
      "position"}]}], "*)"}], "\n", 
   RowBox[{
    RowBox[{"OutputTypeHippos", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
      "\"\<hippos\>\"", ",", "\"\<LinearHippos\>\"", ",", 
       "\"\<LinearHipposD\>\""}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hipposStr", " ", "=", " ", 
     RowBox[{"StringJoin", "[", 
      RowBox[{"DataType", ",", "\"\<_data_\>\"", ",", 
       RowBox[{
       "OutputTypeHippos", "\[LeftDoubleBracket]", "HipPosIndex", 
        "\[RightDoubleBracket]"}], ",", "\"\<.mat\>\""}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
      RowBox[{"choose", " ", "the", " ", "type", " ", "of", " ", "non"}], "-", 
      RowBox[{"stance", " ", "slope", " ", 
       RowBox[{"output", ":", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"1.", " ", "Lineariezed", " ", "non"}], "-", 
         RowBox[{"stance", " ", "slope"}]}]}]}]}], ";", "\n", 
     RowBox[{
      RowBox[{"2.", " ", "non"}], "-", 
      RowBox[{"stance", " ", "slope"}]}], ";", "\n", 
     RowBox[{"3.", " ", "Hip", " ", "angle", " ", 
      RowBox[{"(", 
       RowBox[{"angle", " ", "between", " ", "two", " ", "thighs"}], 
       ")"}]}]}], "*)"}], "\n", 
   RowBox[{
    RowBox[{"OutputTypeSlope", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
      "\"\<nsslope\>\"", ",", "\"\<LinearNSslope\>\"", ",", 
       "\"\<hipAngle\>\""}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"nsslopeStr", " ", "=", " ", 
     RowBox[{"StringJoin", "[", 
      RowBox[{"DataType", ",", "\"\<_data_\>\"", ",", 
       RowBox[{
       "OutputTypeSlope", "\[LeftDoubleBracket]", "NSslopeIndex", 
        "\[RightDoubleBracket]"}], ",", "\"\<.mat\>\""}], "]"}]}], ";"}], 
   "\n", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "choose", " ", "the", " ", "type", " ", "of", " ", "torso", " ", "angle",
       " ", 
      RowBox[{"output", ":", "\[IndentingNewLine]", 
       RowBox[{"1.", " ", "torso_hip", " ", "angle"}]}]}], ";", "\n", 
     RowBox[{"2.", " ", "theta4"}], ";", "\n", 
     RowBox[{
      RowBox[{"3.", " ", "non"}], "-", 
      RowBox[{"stance", " ", "torso", " ", "slope"}]}], ";", "\n", 
     RowBox[{"4.", " ", "linearized", " ", "torso", " ", "slope"}], ";", "\n", 
     RowBox[{"5.", " ", "stance", " ", "torso", " ", "slope"}], ";", "\n", 
     RowBox[{
     "6.", " ", "linearized", " ", "stance", " ", "torso", " ", "slope"}]}], 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"OutputTypeTorso", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
      "\"\<TorsoHipAngle\>\"", ",", "\"\<theta4\>\"", ",", "\"\<nstorso\>\"", 
       ",", "\"\<LinearNStorso\>\"", ",", "\"\<storso\>\"", ",", 
       "\"\<LinearStorso\>\""}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"TorsoStr", " ", "=", " ", 
     RowBox[{"StringJoin", "[", 
      RowBox[{"DataType", ",", "\"\<_data_\>\"", ",", 
       RowBox[{
       "OutputTypeTorso", "\[LeftDoubleBracket]", "TorsoIndex", 
        "\[RightDoubleBracket]"}], ",", "\"\<.mat\>\""}], "]"}]}], ";"}], 
   "\n", 
   RowBox[{
    RowBox[{"skneeStr", " ", "=", " ", 
     RowBox[{"StringJoin", "[", 
      RowBox[{"DataType", ",", "\"\<_data_sknee\>\"", ",", "\"\<.mat\>\""}], 
      "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"nskneeStr", " ", "=", " ", 
     RowBox[{"StringJoin", "[", 
      RowBox[{"DataType", ",", "\"\<_data_nsknee\>\"", ",", "\"\<.mat\>\""}], 
      "]"}]}], ";"}], "\n", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
    "SetDirectory", "[", 
     "\"\</home/shu/workspace/Research/AMBER_Output_Test/data/\>\"", "]"}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSslope", " ", "=", " ", 
     RowBox[{"Import", "[", 
      RowBox[{"nsslopeStr", ",", "\"\<MAT\>\""}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Sknee", " ", "=", " ", 
     RowBox[{"Import", "[", 
      RowBox[{"skneeStr", ",", "\"\<MAT\>\""}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSknee", " ", "=", " ", 
     RowBox[{"Import", "[", 
      RowBox[{"nskneeStr", " ", ",", "\"\<MAT\>\""}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Hippos", " ", "=", " ", 
     RowBox[{"Import", "[", 
      RowBox[{"hipposStr", " ", ",", "\"\<MAT\>\""}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Torso", " ", "=", " ", 
     RowBox[{"Import", "[", 
      RowBox[{"TorsoStr", ",", "\"\<MAT\>\""}], "]"}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5206980829927273`*^9, 3.520698089238552*^9}, {
   3.523978947801125*^9, 3.5239789588067827`*^9}, {3.52397911342336*^9, 
   3.523979114631843*^9}, {3.5243169653528748`*^9, 3.524316973652472*^9}, {
   3.5243319180886717`*^9, 3.52433193045708*^9}, 3.524332052043282*^9, {
   3.5243336818168917`*^9, 3.524333682086383*^9}, {3.527031289451779*^9, 
   3.52703129229949*^9}, {3.527080033657362*^9, 3.52708006144129*^9}, {
   3.527080133566532*^9, 3.52708014024951*^9}, {3.5271009843940268`*^9, 
   3.5271010408501663`*^9}, 3.527250709910329*^9, {3.5272507404573*^9, 
   3.5272507718199787`*^9}, 3.52725082256019*^9, {3.527279016037776*^9, 
   3.5272790190211163`*^9}, 3.5272794890519753`*^9, {3.5272832222939243`*^9, 
   3.527283234403532*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\n", 
  RowBox[{
   RowBox[{
    RowBox[{"A", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"a", "[", 
        RowBox[{"i", ",", "j"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "5"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "5"}], "}"}]}], "]"}]}], ";"}], "\n", 
   "\n", 
   RowBox[{"Tf", " ", "=", " ", 
    RowBox[{"Max", "[", 
     RowBox[{"First", "[", "Torso", "]"}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"(*", "Data", "*)"}], "\n", 
   RowBox[{
    RowBox[{"HipData", " ", "=", " ", 
     RowBox[{"Join", "[", 
      RowBox[{
       RowBox[{"First", "[", "Hippos", "]"}], ",", 
       RowBox[{"Last", "[", "Hippos", "]"}], ",", "2"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"NSslopeData", " ", "=", " ", 
      RowBox[{"Join", "[", 
       RowBox[{
        RowBox[{"First", "[", "NSslope", "]"}], ",", 
        RowBox[{"Last", "[", "NSslope", "]"}], ",", "2"}], "]"}]}], ";"}], 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSslopeData", " ", "=", " ", 
     RowBox[{"Join", "[", 
      RowBox[{
       RowBox[{"First", "[", "NSslope", "]"}], ",", 
       RowBox[{"Last", "[", "NSslope", "]"}], ",", "2"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"SkneeData", " ", "=", " ", 
     RowBox[{"Join", "[", 
      RowBox[{
       RowBox[{"First", "[", "Sknee", "]"}], ",", 
       RowBox[{"Last", "[", "Sknee", "]"}], ",", "2"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSkneeData", " ", "=", " ", 
     RowBox[{"Join", "[", 
      RowBox[{
       RowBox[{"First", "[", "NSknee", "]"}], ",", 
       RowBox[{"Last", "[", "NSknee", "]"}], ",", "2"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"TorsoData", " ", "=", " ", 
     RowBox[{"Join", "[", 
      RowBox[{
       RowBox[{"First", "[", "Torso", "]"}], ",", 
       RowBox[{"Last", "[", "Torso", "]"}], ",", "2"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Tf", " ", "=", " ", 
     RowBox[{"Max", "[", 
      RowBox[{"First", "[", "Torso", "]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", "Parameters", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"HipParameters", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{"a", "[", 
       RowBox[{"1", ",", "1"}], "]"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSslopeParameters", " ", "=", " ", 
     RowBox[{"A", "\[LeftDoubleBracket]", 
      RowBox[{"2", ",", ";;"}], "\[RightDoubleBracket]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"SkneeParameters", " ", "=", " ", 
     RowBox[{"A", "\[LeftDoubleBracket]", 
      RowBox[{"3", ",", ";;"}], "\[RightDoubleBracket]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSkneeParameters", " ", "=", " ", 
     RowBox[{"A", "\[LeftDoubleBracket]", 
      RowBox[{"4", ",", ";;"}], "\[RightDoubleBracket]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"TorsoParameters", " ", "=", " ", 
     RowBox[{"A", "\[LeftDoubleBracket]", 
      RowBox[{"5", ",", ";;"}], "\[RightDoubleBracket]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Parameters", " ", "=", " ", 
     RowBox[{"Join", "[", 
      RowBox[{
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"Join", "[", 
            RowBox[{"HipParameters", ",", "NSslopeParameters"}], "]"}], ",", 
           "SkneeParameters"}], "]"}], ",", "NSkneeParameters"}], "]"}], ",", 
       "TorsoParameters"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", "Models", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"HipModel", " ", "=", " ", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", "t", "}"}], ",", 
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"1", ",", "1"}], "]"}], "*", "t"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSslopeModel", " ", "=", " ", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", "t", "}"}], ",", 
       RowBox[{
        FractionBox[
         RowBox[{
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"2", ",", "1"}], "]"}], " ", 
           RowBox[{"Cos", "[", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"2", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}], "+", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"2", ",", "3"}], "]"}], " ", 
           RowBox[{"Sin", "[", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"2", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}]}], 
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"2", ",", "4"}], "]"}], " ", "*", "t"}], "]"}]], "+", 
        RowBox[{"a", "[", 
         RowBox[{"2", ",", "5"}], "]"}]}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"SkneeModel", " ", "=", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", "t", "}"}], ",", 
       RowBox[{
        FractionBox[
         RowBox[{
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"3", ",", "1"}], "]"}], " ", 
           RowBox[{"Cos", "[", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"3", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}], "+", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"3", ",", "3"}], "]"}], " ", 
           RowBox[{"Sin", "[", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"3", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}]}], 
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"3", ",", "4"}], "]"}], " ", "*", "t"}], "]"}]], "+", 
        RowBox[{"a", "[", 
         RowBox[{"3", ",", "5"}], "]"}]}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSkneeModel", " ", "=", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", "t", "}"}], ",", 
       RowBox[{
        FractionBox[
         RowBox[{
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"4", ",", "1"}], "]"}], " ", 
           RowBox[{"Cos", "[", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"4", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}], "+", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"4", ",", "3"}], "]"}], " ", 
           RowBox[{"Sin", "[", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"4", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}]}], 
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"4", ",", "4"}], "]"}], " ", "*", "t"}], "]"}]], "+", 
        RowBox[{"a", "[", 
         RowBox[{"4", ",", "5"}], "]"}]}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"TorsoModel", " ", "=", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", "t", "}"}], ",", 
       RowBox[{
        FractionBox[
         RowBox[{
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"5", ",", "1"}], "]"}], " ", 
           RowBox[{"Cos", "[", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"5", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}], "+", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"5", ",", "3"}], "]"}], " ", 
           RowBox[{"Sin", "[", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"5", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}]}], 
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"5", ",", "4"}], "]"}], " ", "*", "t"}], "]"}]], "+", 
        RowBox[{"a", "[", 
         RowBox[{"5", ",", "5"}], "]"}]}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"TorsoModel", " ", "=", 
      RowBox[{"Function", "[", 
       RowBox[{
        RowBox[{"{", "t", "}"}], ",", 
        RowBox[{
         FractionBox[
          RowBox[{
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"5", ",", "1"}], "]"}], " ", 
            RowBox[{"Cos", "[", 
             RowBox[{
              RowBox[{"a", "[", 
               RowBox[{"5", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}], "+", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"5", ",", "3"}], "]"}], " ", 
            RowBox[{"Sin", "[", 
             RowBox[{
              RowBox[{"a", "[", 
               RowBox[{"5", ",", "2"}], "]"}], " ", "*", "t"}], "]"}]}]}], 
          RowBox[{"Exp", "[", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"5", ",", "4"}], "]"}], " ", "*", "t"}], "]"}]], "+", 
         RowBox[{"a", "[", 
          RowBox[{"5", ",", "5"}], "]"}]}]}], "]"}]}], ";"}], "*)"}], "\n", 
   RowBox[{
    RowBox[{"TorsoModel2", " ", "=", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", "t", "}"}], ",", 
       RowBox[{"a", "[", 
        RowBox[{"5", ",", "5"}], "]"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", "Fits", "*)"}], "\[IndentingNewLine]", "\n", 
   RowBox[{
    RowBox[{"HipFit", " ", "=", " ", 
     RowBox[{"FindFit", "[", 
      RowBox[{"HipData", ",", 
       RowBox[{"HipModel", "[", "t", "]"}], ",", "HipParameters", ",", "t"}], 
      "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSslopeFit", " ", "=", " ", 
     RowBox[{"FindFit", "[", 
      RowBox[{"NSslopeData", ",", 
       RowBox[{"NSslopeModel", "[", "t", "]"}], ",", "NSslopeParameters", ",",
        "t"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"SkneeFit", " ", "=", " ", 
     RowBox[{"FindFit", "[", 
      RowBox[{"SkneeData", ",", 
       RowBox[{"SkneeModel", "[", "t", "]"}], ",", "SkneeParameters", ",", 
       "t"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NSkneeFit", " ", "=", " ", 
     RowBox[{"FindFit", "[", 
      RowBox[{"NSkneeData", ",", 
       RowBox[{"NSkneeModel", "[", "t", "]"}], ",", "NSkneeParameters", ",", 
       "t"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"TorsoFit", " ", "=", " ", 
     RowBox[{"FindFit", "[", 
      RowBox[{"TorsoData", ",", 
       RowBox[{"TorsoModel", "[", "t", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"a", "[", 
          RowBox[{"5", ",", "1"}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"5", ",", "2"}], "]"}], ",", "10"}], "}"}], ",", 
         RowBox[{"a", "[", 
          RowBox[{"5", ",", "3"}], "]"}], ",", 
         RowBox[{"a", "[", 
          RowBox[{"5", ",", "4"}], "]"}], ",", 
         RowBox[{"a", "[", 
          RowBox[{"5", ",", "5"}], "]"}]}], "}"}], ",", "t"}], "]"}]}], ";"}],
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", "HumanParameters", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"HumanParameters", " ", "=", " ", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"Parameters", "/.", "HipFit"}], "/.", "NSslopeFit"}], "/.", 
        "SkneeFit"}], "/.", "NSkneeFit"}], "/.", "TorsoFit"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"HumanFit", " ", "=", " ", 
     RowBox[{"Join", "[", 
      RowBox[{
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"Join", "[", 
            RowBox[{"HipFit", ",", "NSslopeFit"}], "]"}], ",", "SkneeFit"}], 
          "]"}], ",", "NSkneeFit"}], "]"}], ",", "TorsoFit"}], "]"}]}], ";"}],
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", "Plots", "*)"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"ListPlot", "[", "HipData", "]"}], ",", 
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"HipModel", "[", "t", "]"}], "/.", "HumanFit"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "Tf"}], "}"}]}], "]"}]}], "]"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"ListPlot", "[", "NSslopeData", "]"}], ",", 
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"NSslopeModel", "[", "t", "]"}], "/.", "HumanFit"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "Tf"}], "}"}]}], "]"}]}], "]"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"ListPlot", "[", "SkneeData", "]"}], ",", 
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"SkneeModel", "[", "t", "]"}], "/.", "HumanFit"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "Tf"}], "}"}]}], "]"}]}], "]"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"ListPlot", "[", "NSkneeData", "]"}], ",", 
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"NSkneeModel", "[", "t", "]"}], "/.", "HumanFit"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "Tf"}], "}"}]}], "]"}]}], "]"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"ListPlot", "[", "TorsoData", "]"}], ",", 
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"TorsoModel", "[", "t", "]"}], "/.", "TorsoFit"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "Tf"}], "}"}]}], "]"}]}], "]"}], "\n", 
   RowBox[{
    RowBox[{"matrixName", "=", 
     RowBox[{"\"\<Afit\>\"", "<>", 
      RowBox[{"ToString", "[", "LineNumber", "]"}]}]}], ";"}], "\n", 
   RowBox[{"Aout", " ", "=", " ", 
    RowBox[{
     RowBox[{"A", "/.", "HumanFit"}], "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"1", ",", "2"}], "]"}], "\[Rule]", " ", "0"}], ",", 
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"1", ",", "3"}], "]"}], "\[Rule]", " ", "0"}], ",", 
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"1", ",", "4"}], "]"}], "\[Rule]", " ", "0"}], ",", 
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"1", ",", "5"}], "]"}], "\[Rule]", " ", "0"}], ",", 
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"2", ",", "5"}], "]"}], "\[Rule]", " ", "0"}], ",", 
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"5", ",", "3"}], "]"}], "\[Rule]", " ", "0"}], ",", 
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"5", ",", "4"}], "]"}], "\[Rule]", " ", "0"}], ",", 
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"5", ",", "5"}], "]"}], "\[Rule]", " ", "0"}]}], "}"}]}]}], 
   "\n", "\n", 
   RowBox[{
    RowBox[{
    "SetDirectory", "[", 
     "\"\</home/shu/workspace/Research/AMBER_Output_Test/data/\>\"", "]"}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"matrixName", "<>", "\"\<.mat\>\""}], ",", 
       RowBox[{"\"\<a0\>\"", "\[Rule]", "Aout"}], ",", 
       "\"\<LabeledData\>\""}], "]"}], ";"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Export", "[", 
     RowBox[{"\"\<Afit.mat\>\"", ",", 
      RowBox[{"\"\<a0\>\"", "\[Rule]", "Aout"}], ",", "\"\<LabeledData\>\""}],
      "]"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5206980829927273`*^9, 3.520698089238552*^9}, {
   3.523978947801125*^9, 3.5239789588067827`*^9}, {3.52397911342336*^9, 
   3.523979114631843*^9}, {3.5243169653528748`*^9, 3.524316973652472*^9}, {
   3.5243319180886717`*^9, 3.52433193045708*^9}, 3.524332052043282*^9, {
   3.5243336818168917`*^9, 3.524333682086383*^9}, {3.527031289451779*^9, 
   3.52703129229949*^9}, {3.527080033657362*^9, 3.52708006144129*^9}, {
   3.527080133566532*^9, 3.52708014024951*^9}, {3.5271009843940268`*^9, 
   3.5271010408501663`*^9}, 3.527250709910329*^9, {3.5272507404573*^9, 
   3.5272507718199787`*^9}, 3.52725082256019*^9, {3.527279016037776*^9, 
   3.5272790190211163`*^9}, 3.5272794890519753`*^9, 3.5272832222939243`*^9}],

Cell[BoxData[
 FormBox["0.5867562797946462`", TraditionalForm]], "Output",
 CellChangeTimes->{3.529855415879889*^9}],

Cell[BoxData[
 FormBox[
  GraphicsBox[{
    {Hue[0.67, 0.6, 0.6], PointBox[CompressedData["
1:eJwVkms4FIgChodFbbNEQtZDy3Kao7ZtzlNsqf3ayCXVbFSzSS2xbpVCJ8dW
tlApFBkqtCzLg1yaGddcMhODGcPcx22YcSmWEJUwy2l/vM/3//1eyzMXPH7R
JBAI/p/4Z39ujmy9vikIeuyr5CpTN7iIH6flmMfg3FIapdUhDJGVxxZ8yPfw
gOV0gD71G0S8+LJkqzTI3lLCxxZuYaI+QHlxSyY65+aDbEvvQYNqdYeUmYNr
5KgvuEY0JFEMcjfk5IFIefJKZfQIX7oS6lbnF8DvuDhmvWsW8vdOiz8WFcHE
M5cV3JAD8s6BybGyEoS+qdFYsM1DI5mv3cMsB22hIKb4Xj6sjU7f9/em45lt
YcGScyFCUqSZ/QMMcE9fm/EfKQZd91DhUd8KMH6vVCe5lGL+TnNFx3Albltv
dt8xX47dOrtZTgHVsPnLasHjOB0xMUx+/VgN7PsL6DeKGGhftu3dfvY5mFYK
RrFmBfSu5L4qeVOHzJc60+PUShydN52zvtiAnfHCrrCnVciISFnJmm0EYQuZ
8G6hGjORWrfu7moCKcjwcfHJWhja7U3kLDVBeePo62d1z2H37soDzQYW/i7R
Oyw2qscJRvWj76PZMNoe/HxrWAOuXpz7/Ve8BDHX0USjrhHZW7/NryI04z+9
XduuTr0AezLk6SyrGQZmzKp6WhNGiwvoW2NbQCyT3OnaxsKq4KHqECcOfvjV
gvqYw4LtJovGAu1WKM7OXrrsx8bh0RPNQ5xW+FU+6i1dYiMsL41rEd+GerWw
ojnhJWi+QoGXWztMvTV01lk3o3qjrjx9DRd/+nxpoMFoRo/CVSHicfFN9WJn
qEML1Jlxw3pJPIg+ezbrzm/BRq+m8QOHO3B5+QO19SAH+zaop2+t5cM82uml
I4+DAJn9B7aAj/AeJ/3dJ1txhxahXk7pxNqP4qTbylaUeJRrOnh2wbmhQuHq
1YbA+P/GOawRwD8i1MFA1IaowNTbPSECjCg0V094tuOuM/1uJE+APcEWy/n8
dmTZdCUZbRFi5fCoo+d+Lkq13iQzE4XI2Be72b6WixfDa2hH3ggxEcxIMrPl
QcAmPZw+JEIhy3370wweVH84ZySViSD8XBpzX78Dc9f9n2xeK0Zi+Tcd393o
gLZPTE77BTEqseW0bKIDxsjJCxSI0RiqU5zqzQfpH91kCY75tjWtbuZj5999
RXkpElgF3+BYkjrh3r9Q8sOsBLxV0rhHqZ3wrjN5NughhdQ424W32InzGTuY
15hS0DZc+0nrTBeiozyrzNbLoN933/ICtwvJP4XV1l6SYU3UiszMWIBc+/v1
VKkMThtZXvHOAjCNS1+83yGHjR3h9VexArS857JT0+WoLwlKGqgRQCYZayHP
y0GZD0/h/iXAGFOnvYvajdoi6tFAcyEWH1h3nK/phpqY6nvzmBDE8H1dRNMe
WHjGuWbEC2F+xEdUFNWDlIeVn6U1CLF1W7TUpbcHs+MJ2Tofhdi7Nqt7dFcv
KI6mIfZkEX6cqu2LzezFEbMPX0yeE+EMXz5gqe7FMaPzN8MLRYgoea964d2H
P7qjFCXDItxMMBw91dAH81Nages3ivEwhDy2ZN4PtxHeIO+kGIVulInH0f0g
hQZ9m50uxnPS+Sn7wX6Irx8fPyAWf/Ka8FYKBbxf2yy9XieB4lXhu4gcBZJn
356opUgw1cKZN9AYQKPDQgo9UQJC/shiue8AqKxBlRtfgnVxmsuH2ANIyO88
JNWV4mu/rzQmrQbBuxh42p4ixfZ932vdjR1EdXzN/tZkKfZbeq8ijQyCnq7O
FIul2EzYo+1npUSk/nhRnKEMFdM8HYaHEhO+pJmzXjLsGfRarRGrBGPTlOxt
lgyczvHPf2QqP/Wq9HbplYHS+D9i9rASfVm6lKwNcshLV+lOGapwbtJDK/WU
HD5P0vX2OKmAf9uhNkuO8UQb/cRLKqgv15HYSjnCr1YY9P2pAvHhdcsiUjcW
zzoa2kpVoB1hlH0X1I3Yk6L1UdpDUKWZaVMLu0F09zVu2zGEgx9WdA/OdYO2
a8bEJGAIluH/0vrargfmtr+ZBqQPoa++v2Xnp98LTPXMKjlDKJvZdGWlpQf/
B47mMMM=
      "]]}, {{}, {}, 
     {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwVxXk01AkAB/ChXk9JKFcKScmq6HTk+DJ5SVJ2IqSknBstkWZfy7aaemTj
zTDGUT1K3kiOF6EojSPNyDFjZpj5YfCzzmezrrYUdvePz/sYX46ihShTKBSP
//z/zc1lXIMaN8dvguq+3ZuojnESqeY7zi6o5eXKUtjbsPlubsEfHGvkj1Ar
Qtn7ENEfP3WN44oJFd08L7YT2Gs/aYRwfJCmn5jpzvZEUZ5pQQAnDIEm1Rqu
7EBUTFIj/Tl08LmeHS7saFhriIR+nCQsTRiUHmX/jrio6D4PdjZGj7cWNh1j
gf6T2YOj1CcoiOPOTxfngCZySQvdWoSoNIug+Nl8RFxU+8SIK8O8U4lGQ3gh
WrUv8XT2V2BkYHGpNf0ZltI0eeOhVXhzYcpSW6cMTiYy8m3qK1xZXygypb5A
SreWkd9vdTAPtH++80wl5h5yWdfK6nHRZCe3mVmFA86yDWYSHgTnP1OPFdUg
piFsNNutEV0vE1Ycq17D2X27H1/RhFPxZxU9n+rA/jNLRvi+B6VW98IZ5XrQ
k2zmns+2IOOg9i1rYx6GrIPzlu/wcSPn6wDFvgGd/5zwPfRdgNgjXcy6o43w
aDhfYpbxEbszxLSWS004nfuDGkW9HUumguDGyGYMxhuv5RR34F7MWNKGtPeI
t6J9PGskRHB7oJ1OcQtOWF3Tk/QKcS6ow1G98gMGa3dW7isSQfVF8BnfLj5s
PlhLlkO70G07XZOuEMDrnHq92nYxZkSSkuSRVtjExC98mxKjxXi4okO5DQpN
10r/YgmqS0ZU1qi1w92KW59Ll+LitEZ06t4OxA+1pes7d8PXOStMeLwTUVqe
f6kudcNqLNnntpkQ5m3hmX38Hjy6RzCaa4Tgqg40GSbJwHIz9b3lIMK4SL/u
sbMcFqy8+SyeCH9vnVpZtZ4A91BpjuxUF2qn9T/28AmULrwcZYm7cNJ1YDgl
sxfMjT/axwaJkW+qxCj36oOSx5ht+4QYr+6PT2pp9ENf/fnKnasSZNND8v16
+6Gxd1VH8HcJbCfRp56jQHq4t8vbRCmeDJWHFIYOYHbFrslhYzfGy2d7HfQG
sWjouW5rfjcGD9jzL90YBNM/4GvMjh7Y/eLpH9c3COV5ysxMZQ/oiRkvGuyG
oFT8iDlhK4P2KdtIi5whFI2FfXF8K4Nu/7P2lOUh7ElsCX/jLseipdus0QUS
DjZGiRVCORj3u461BJCo7V97aLVYjvUT/g8iA0lYM+ZHvaVyGD752eV1EIn9
HYKTX+VyULXYHNoVEjuCr+s5DcuR8kVx5C6dhCpTUN6+IMeWhuuMSRYJYixW
Ma5P4KnBcg8zg4RPagDriAEBi5tJe6wzSUgPuLncNyLgdDBXysgm0ZlgWGy5
g0Bw4btdBnkkGjcJ4m7sJVByb13b6RISz2CottqJwOHRjG2fS0mYjqjwvKkE
6qkGcQ/LSRSkzMVwXQgIv1kaTlSQeCTly064EZi/6h19+zUJZkTsUxaNQEKr
otmsjoSaZoDPsBeBNbvCN3e+IZFSfXzdYR8CeoM3G7fwSDCUDKNl/gQe26/W
bWwgQeGqmJgHEDDPSY0IbyKRcHJO+msggcoFHd6G9yQWZ/qT2y8TsKfla1W1
kKBn8e2MQgj8C9J4bIQ=
       "]]}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->True,
   PlotRange->Automatic,
   PlotRangeClipping->True], TraditionalForm]], "Output",
 CellChangeTimes->{3.529855415922308*^9}],

Cell[BoxData[
 FormBox[
  GraphicsBox[{
    {Hue[0.67, 0.6, 0.6], PointBox[CompressedData["
1:eJwVkms41QcAh5FDEUIyJOvQ0xk9yiKr6GcNW7eFNDK1o4uaZrUkdFGjcglh
pQvqLCOaRI5b4o9c4jgc54Jzca5/hXKJRFLWPvye98v76X1+y/cf8zmkoaam
dvDz/ucvjREtF1YegW1hqtg6RYjvebeuMyxj4HimkjOzXYSIst3v6Q4p2G3t
XW47IgKXFV+USr0OL0pdOeNPMV49DZYfX5WJ3/byzanmEqj7URNomQxEb13H
tm+TIHmn4b0vGDmoXSZySD3UB/Mf1Krn5+ZBM3s0aQ1Fily3Ud50QQEYli53
zldK4bBe+nqgqBAmnm9tJ/bJUOvApghLH4F1QO+3yDkZbEz2XT0YWILAh/p3
S/bIEZImyJRIH+PU+UUaD4rkKNHbke8bxISvjfmYcE6OqYRGZruqDBcTxD51
vgq4aLnUuwdXgLlw7wL6fQViYkrZTwcq4XiFzhVNKtD6yVbkePQJDAq9giM8
lNA/c+9F4XA1qPS4Uz0pSvhOmU3YHK/BkcVCmhlXidthaXNZ47UIsR9NLDBS
YSxC83LihjqwTK+0Be5SwXidW1Lzhzo8T718uCJFhXVvz6Rr1NSj9ZubIr9n
Kux5XHFzU3QDokteyXynVDh7fOLOaTyD61zLiSEaibv2q3PL1RohV9DLVnuR
aHgd8u94fSMOKuKoTlEk+h/kldjHNsHSR00UkEFC+1dlRYh7M2ZX31/ELydh
u3JZbR6lBZ4rztQxuSR+7N/TqGxuQfc78/7gQRJ/5FxvWxb/HEusTAz2TZG4
FtTFCdjSivq19v/0zpCosNLrydBpw+pD39BejpMQ9v3Qx2W1wUV/nnu0gsRs
5kWVfjILGz3cqxmdJKwC6ga3/tgOD+fshBcPSGz+Ynb0sgEbkshOfn4MieBu
53cNHPbnroF0HzcSCdfCZj+ldeDPeeuYTW9VKPR5pLFxVyeSBdq1Ov+qcDg+
/OJGHQ6+pZ2ODw5QIerwX3HCEA46qvXWSz4qkehZkhjB4mDhvXMzWTeVyFrR
mWyyqgtHyzq+3PW1Eg81h1NLk7qwdxrpC2sVIFQ617yHu7Bfok7f6aEAp4F2
Y3QHF88EdtVpzXIo/va8nVzEBUlxyvR2lmPiwsFsOwMeGpPHeSO+MlDoMYzW
Yzz4uLiltLpKsQSMnMMcHvxCjMI7NPtA+z+3Ax+O83RrLHrFWP9RXJCTxgf/
w0+DxQwRtkneF347zkcg1SjBP0iIwGrTYpmPAGf/trN+T+tF6G2n0nOlAoQX
b2WaSbsRHbWr3GJxN3aP1rTePidAqv8fVVUnuxEzP2vuqRMf95yvPvUTdOPq
YNTalWNclC55SEw69eDNiLXn0OUuNE22NfyV0YMjM+yT/cs56OYPNDlM9WCp
Yd7PxVlsDJRqtXb69cJmf/40fbwNM+k27aGVvSjTCR1MmmmB7onNnbpmQlQ8
eFLit/7zz7zp3IIoIerMto/KLzXAfk204HuREF+9LIl7F07AzSCrt3+DCK+U
pLHVpyp4jVSJYzNFuKDlaHi0ion97B7p8lkRit/rFam0ihBWOKkgAsWoq7Dg
rEnJwaUrxv17a8TIoF7YIdqUjhshDgMfLCVw4xqlWi77lcjfsvPVrWgJ2nM+
Bp3zvkE8oYWOOMsk2OLXd3r4p1yCpX3ljQB9kKomv7IeeUT0vch/G8boQ/iC
U8c6h5nESFPzlKG6FHODRnHvnlQRarnkzKMgKVYFrb3DaqwljC5qfNrRIMX0
icbsoaZ6wvrAl+qvqTKETIdXFxs2Eo6bN2kmxspgJZiJXvV7M+GxPFCbRsqw
3T/szb6q54SdmivlAFWOUX3/yYmXbQRzlKX12EeOagn1O7uP7YSrLGC+eqwc
mvWF8ZcmOojmjsEFXqVymBKRJnu8OMTO2kjduyo5Itv8K+9bdBE9D7X1RowV
qAwOiovmdxH07Ax9V3cFhFTa2PYjXGIwacWipJMKTKzJ21YzxCVOnGUaiv9R
oIo5ns/z5REzR78zthUoENqz1H8zg0fE/sxdHEVRwnJp6G62mEfobgta8txJ
CcWtdEP7aR5xbcOYqWmwEuFS9aEXanzC0va8WXDGZz/nrDHlDY/IM9O3KGtW
IjA7K8m7gUf8B3JVRW8=
      "]]}, {{}, {}, 
     {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwV1nk8VG8XAHBaVMgWUdmTUpSENhwhochPskTIXpIlIYQiIVuFJsmekBAp
S1zrzNjHXHu2uXOTVFpIUvQ+718+34/ruuc8Z3lkHL3MXFZxcHAEcHJw/P9n
0JaSZxJvDLW+5B40jhEZgmt9/YJY6k5waKunT20bgi2303Lvph4EbYuqQD3Z
IfAYC/nsk3oCij2ifTn2DUHyhlkBl1RLKK7SqAwyGIKCTPlcu1Q3iE67UREQ
PATlMzqXbVIDgOom6x3PGoKDAr0M69Q78LknznbhxTBc8/IeNU6mwCpL6XEH
93cQcHHXY12dHHg8x3lXknMMzHr1ElzFCyA/4UD6dMw4eNhvnI24VgInigzT
7DZNQrvIhYbN+8uh8YrRdE3HJCwnCDZMu1YCXedwhXI0C7S3DxF18VXQ3kFp
9D5AQOyAsJR1aC28ZjvPTbQQMJf+7J5PST2sqbDJGDNlQ3TDII91TT0UlEh2
NJmxQZJcF6VNrYcXX4WG883ZYLTHPYh/vB6Uit6lelixIad6l9OLjRhkZXAn
fLBng+lAoerHyxiUFm/lLbzChoF1Ax5lIg0QXOhtcfouG1SODfHt6msA15jd
fwQb2XAgTJ25PNwA6mryEnlNbFCtS07BJxqA//IuIbUWNqgf/k887FMD+PE3
Xj1DY8OR/e0KA6sagSOqi/tWFxt0ZN7qRe5vBMq3NWr5w2ww48wKmkxoBAIL
f7/uOxvOwIrm6+RG4E18W3b1BxvMb9hyxqU1wpugo7rjc2yw+C0WfTC/EZo3
339busCGc9/vpSTUNULanpQF3b9scGRFlGl8boSvHF8P7l9Hgm+j2xTFsAlo
G0/4i4uT8Oi5St6X000wqVv3M1iChIaU5Qs6Fk1gdu3Nr2FJEvgvPRj75NgE
Yl7Xbt6TIaFYCOvTCmkCX83ezFl5EqacNjeTJU3QJO2R77CfhHNcrZnKws3A
cb1zUPo4CTe/JZ6/vbUZJEPW957WJ6Fg5Ny2Eelm+HuBD26cIOFXybfUCKVm
MF/XFYAbkpBsKZ4woN8MRm/9TD1NSOguuBoScr0Zhqa5dCwsSTh2UtaaPt4M
jTlKxktuJOTLSOeWkc1w5bHEwSl3EngXJb5QZprBib+8jnGRhIG8LTcvLjTD
t8WbmdkeJFxeESjk5msBDYvhYhUvElJf/ls8pdkCLxqanwn6k/DnzrKOmm4L
SPXzqH9AdrD7Eydh2AIc0d9sagNI2MPzS3rWvAXmdM942l9H+XKeNUi83AJx
2hNPUkNImBEdo/Q+boGBiMzvDbdIOD07QlRnt4Drky99YREkvGoZUsx51gK8
GYtGmpEkhPn0NfhWtID5stBsxW0ShDs6pjd1tMDfYc71SdEkQGjNIYulFvgQ
kGKwOgHlh3w4NGLVCvkbKXP4QxKk3VZJ6Ni3QkN50UtvCgnPP16+UOjSCnmN
z97wPCKhafbYTIBvK9gv/XykkUaiPHz6IxzXClkrnCUx6SQYb9SWOo21QumH
1FZGFglDiUVOr1tbYUOE1tB/2SQ4C4kUSHSi3+uDUi9yiOiM8uehVhDT/Zbe
nkNCkUyyTsyPVvB1eBRTkEcCl/q0S/MOKrwmxI0FCki4X2VWtFuRCpKcfOev
IUscrZu9p0KFN3vWFA0hq2rfD3AAKthLVvx+VEiCo5FGzLIVFawiBN3WPSdh
tjO/y9meCpmWUhZ2yNdPCwp1ulBhKi4zrAI5yXwqLc2XCmZ+5eFWxSTU2yUV
H4yjwmTiaZX4FyQYsZa+ZdynQg23+Od+5H4nFzWuR1TgERgaEy8h4bP7kfq+
p1QY4iHcniIHfspbpVFMhbST9/9NIa++wn8it5wKCQULXfKlJGy9Svb4YFQ4
LpC7MRs5/6eJ8HArFdzdToe+Q94fWG2l3UmFEr9AOeEyEgxuJLD4hqnwQZBX
IBw54M6huec/qKCuOGij85KEQ+ofNX8uUuHwc8ZTT+Ql8lG01j8qqKk/OvAQ
OVTnj0QvDw2CxIO3vEfW/v7cfasQDXZPcflsKCeBM8u2wkmMBp0UfR5F5Mjl
OoMFORocUmvjuoysX3zlAeyhwU7/cLcY5PU2UuPR+2lwy6Rvw1Pk9g2MXcyD
6PmIsC/1yHerwq9u06JB5WZNoUHkU2776531aLBoQA/4gsy3mVhfYkSDrIej
EqsqSGC03D/zy5QG6+rV1ogg37+qm6FtSQOzysTd8sjmsvPTMedpcPtVfpIa
skhv3gHciQYU5YPqusgDYWdDxS/RIH50eetpZMperjYXbxqEFL7RskY+N/Z6
U6k/DcqpMtkXkMXj3OwWQ2hQlf37mDvy+BGxwmMRNPCf+ivtiZz5kT4XG0MD
tbCPmt7IDpTrWn2JNGhOvkHxQZY9sTtGIhXlMyFA+f8mf47gruk0sBVKXeuF
nJ93V7IshwbyaymiHsjuZzQu/i6ggXWshqsLssKqLxU6pej5NbrfzyPPlD1Z
uVtJg1hhjypz5GJ7E8P+WhpoP/arNUS+wvfvgWQTDRbeSC1pIO+rKx13o9PA
vE3Rby/yNw8HhZfdNNDytFWQRC7fKui31EeDH+FeW3iR/doa63Xf0UBAZqfe
Isq/eqDvhngWDV4V6hUSyNX9eIbULA3GZatkypGDIiM/us/ToGTDfdWHyBoH
1FTLl2jQ1DocHoTckJTapreODjvd61qOIt+CE8IJG+lA0Z2o3YKsO/vLbnAT
HcrGgz79RPVGPWk9f1GaDkf+Kc0VIEcvbYAKeTpwBne1hyEbFdbE/FWkg1bu
wrA5cheXuFTiYTooFy0++Y3qPaGy8+IQ0KHzlo1ZG7Kp841XMvp0GPpjo/kQ
ua9x3PCVGR1KYk+9VkJO9U5MXraig2kWXXcB9ZOllPaEvj0dVC2WuOqQR0Ky
/YY96GDk+UjpODLroHPmSiQdIoWTfgLq10nmCf93sXSYv8zzZQX1+6TnHuOq
JDp0k67mb//v3B9LPul0qOnm9lVBHue/afG+gg5CZQokD5of40XOSk3VdLjW
RTtbj+bL+HGDNZkYHapMSTUv5LFg/nKrDjoUO3zS7iwiYfTDk42dBB2iIsNv
XEXza/TWTbJgmg6H/rO32/J/S7jU3p6lg/Mpt+I6NO/enVG8CEso3tmTiZzI
I1hNS7lgG0j0ER7++Wj+PhwMfqTdBueWW2yFc5FVas389dtgyVZwOQvN26Gu
DIUzp9rAzn3VekXkwdWug7xWbcB57tB5TTSfB7zmVcK92oCXRuMwzkT5PSE4
45rRBs2p9TzSj5GJ+QbdvDaQ8yzVe4D2Qd+NoYfSRW2QnDP3fg0yXpF5fKSy
DW7NTDuw0T5hSu3NNu5qA997larxqST0/DKyUv3bBj+e/lJLvkfCvNZjhaRV
7eD1m6ifSELzMerT0qf17SAmW+uigOwqHPckV6Qd9B8m8L5B+2tFuYsQ2tcO
u+LO1dTeJWHvpdOe3y+0w6U/R/wuof1n/jJTy9i9HcajXdseov0YtPiVv/BK
O1hXCR1pRvuTGpVU7hDcDmvW2HoKo/1ql9v7qye5HfpedYhmhaJ6Gz1zq4TW
DqMvxs/bon39anveGe7udtj9u+ycB9rnI5fm5Vz70Pe9YXkHXiNB/ncyTYLV
Dh/UJ6fjrqL9ITLAG7/UDgqeX1ZnovvBrInVQw+lDtDTlul0Q/cN40bb4l0P
OkDZTYMdexb16/zEBrFHHfByu/VqH3MSNu9ycluX2QGM9nuBZ8+g+km4KDtV
1AFV3/kct/yH5oGtPyW3qQP4JEz2xZ9C3/srMVLyRwckqp7Kk9NF9x2lZlvh
/zphXLI+UWgfCSoX9GtWW3aCbM6IarMS6o/kNtE52064tEVPyVcRxf+nh9nr
3glad3HRDgX0vfRRg8TwTjCSbWO6y6F57bygyl3WCU6Ti4NntqD7SJrCRg7+
LrhSvPHTB06Uj+tbXK6IdIFOQsoZHQ6076031I1u6wLHO2+b01bYECj28XL1
ri4wvTxUc/wPui+nFnT66nTB1m1fV0fOs2H+/s6499e6wGjT+3/ZU2yg3N3B
0zHaBWs3l+Y4tbFhMkRmQ2pRN+TYY/SVaDZkRVu6TZV1w6mgX4rvo9hwITm+
Vf1NN0QuJBe3RbKBeP775mBzN6y+qLMSF84G9kjvkthYNzSqGn75G8CGqUO3
Pj/m7wH5J1Lct1zZ8Hme6Mn264FbUVYMGR02hKibdVhIMWCZ/+aOujkCgsuj
GgxkGZCrW2Ni+Z2AoL21lUd2MEB4Zn7w6ywB1+XlsiT3MKAzc4pbfIYA/80L
flPqDEjtFVxwmCTAe+GRpL8JA7Qm1y6ndxLg9Jrlk3KDAZs/qb0ezyXASN1H
rO8dA7ANWTNJxgRM1uyoUC7ohYnCzMzqEBYcoh3sW3Flgt53inRt5iSYn+Ov
3yiLg3S2rXae/QQc8g35+eczDpRkobUy8mMwLniiwqaoDwQGNbUlp0fgpPqz
+rSAfnjgeXy/kOkQhLA67289NgBh6x0/dGf0g5ew6Ree5QFQKXPfQPnHhN2d
7imj9EGwMuE+3sLHgGc8E82Sd4bg7xOpDfw32mG6d2tt9rFhuHTzaPscuwW+
iX/+t5p3BIQKfk9dSaiDmq9bOwbpIzB6DZe/uq8cTp2YYMemvIPwnp96Rpso
kCXPGVFqPgraWdyh8cbJWFXc9IywwBhQnuw0mN9XilECXLKs342B1UuvP4eC
arDDMzDK/2gcjFwb9yWtbsJyWKUuT10nIMMyJMqgg4pNl/54pyk2CdHbw2Qj
uDqwSRUN+gX/SahixGVX3OzGFiKiNy8FTAKlbtT5bWI3xtPf53z/+iSccvah
1T7pxtQDPDmbQyZh/eN8m4zqbiyuNuOIXMQkqBopWP/41o0d1ON8MZU4CSHP
9lxlnu/BEs5S73kUTkJkqNqgohgDOxpoanNtdBLeFtwRGH/HwBwiC+N+jE2C
ntGWqeMsBhaZtKrea2ISCgB0CqcYWHfBK+lLxCQIyVRHO35nYI7DYlN205PA
xZWyM5OrF4s9wvI2mJ+E15SvAaDciw3/9YnaxssC+cVUXCi8F1tZ3/GGspEF
pxyl8OnIXkxWRO6jCD8LptfmK9XE9GKXlQZOCgixIPT3uVmTB70Yh90RwbVi
LOjuOnxdO78X24WtSp/dzoLwk2LxOh29WMDNBy8bj7JA8UxbTLUgEysSC3U+
qMkCizlbXrYIExstdRd9oYX+ftOhb+u2MjHtcc0blGMsUF27PKAry8S4j04b
eJ1ggWOry78IFSaWMa8xKX6GBffF3MvlzJgY1W2KL/ASC1ScPb0n4pnY4gqj
8YsHCwr6Tfp87zGxPam1fk6eLPg7Iu/HmcLEklqSRky8WeheuN9JIJ2JnZM5
mr/DnwVDVwUn5guZ2OxIoiZ+kwUh4Q/NBluYmIjJ4ct7H7Fgd1er+9ICEzMT
wMaa01ggyysefWCJiSUw9U5bp7MgVvQFy3WZiXFZmqpEZrIg1Ytb6+1qHFuw
d1scfsoC8vmzazICODbgkxJ5u5wFVxpiy/t24piQ6raFra9YoCNph5ftxjGT
hSy3skoWaBjOit1RwjFqULHRuyoWpBH/eGUP4NjriGaB/RgLbjH2+glq4Vhq
yvf00Q4WjFv5bes3wzGmZQCfbxcLLr1O/GR2Fsf4ti6HretB8UJxR7sljkVl
rHdSYbIg+l/f6wJbHPN/JqUQPcQCP8koEQlXHLOoNq488B6dl/eq0akAHLsf
jMu3TbHgnO8NkbVBONataU2xm2bB+pi1JyVCcEy/yTk45hMLOsttnh4NxzH1
juBjE99ZsEtU7ue+aBwTHSvsil1hQdTczNzPVBxzPKchu5aDAEnBzy0pFBwr
Gez2D+MkQP6Aj4FyGnofc07Kbw0ByqZHnM4+wTE/mqbveW4CcmxMuuxzcaxB
j0Ed5CEgvOK7aH8ejvE0OW4z20iA2FjCxuP5OJbz9k6LvgABwt7LP4ULcYzx
snez8mYC3p6vz/YrwTHxfc4eRaIENBVo7cZKccyteAGT20LA7noLizUvcWwl
f9vFLeIEBBoWpoRW4JhiukvtKlkC0mIDH6pU4Vh01G/7fiX0/P2x8YcNONa3
Ou7V6X0ElLAHVOMacUzqpuSGdmUU315J6eAmdF7BuuX1BwjIaNxaYdCCY6R3
/JqCwwTIdcWV59BwTHlWylr2KAFRBuVpHnQcC/Yof5GuQUBNl5bgvjZUH66D
FveAAF9yR/7Tdhw7ZiNTGHScgNQ7YZb/deFY3FDF8rw+ATOrFB/PIg+d1Tfz
MiAgNOiWa1Q3jnmbevxxOknAyPl3SU97cOxt1/Lp8VMErPdRn1Vk4Ni6k0l5
ViYoPz2fe0qQM45XGhv/R4A7j5xETi+OzTSdyKGZEbA0serOJiaOqWmPLBwz
J4D6ddEnDDm87vLJt2cJOMVSH3uP3HnkX6a6JYqPa5amj+OYWNW9+TIrdJ42
pvtzkJ3U5Az3nEP5j04SW0Re2mf4Q+o8AXYXA88/6MOxiDimPtWOgLjugSuT
yLwfbR5fdiCATPqcId+PYynH338VckR7jMj67IIsmXNFr9qJAJ3VoxYZyM9W
flHsXQiw0Lwz2oOsbHPzy1o3Ai59jAtZQq5+w61T7E4Ar+eQmuQAjukIJ6ea
XUL5PWqx7ihyh7fEp0UPAvr6576aIpt35UOmJ4q3I3PWHnlMQTn5uBc6r1rD
Ve7IrlHV05+8CXA2IxQuIn8ldDTv+6J6NLZxdkIOhM57h/zQ/wt8VmaBzJl+
dmr8GgHdpTUbdZFjF8eP3A4g4BwrJmgX8qaz7ol7rhNwmZtngQs5/eV3dm8Q
iv+wWtg4+v4dfMGHAkNQvyRxiJQil1xaEy8ZiuIF19eByAdp8ayWMPT9dx0c
jyI3bBdV97iJ4uucEf2F8mkYnhUrGIHed4lz4Dkyc1Rh4k0kAXxk9mNrZJvD
FQfsoghwoLe7cSCTKRrRa6JRPb8NPpyFzsvzR+toUQwBQwGV/IeRF0xO7//v
LgEU+vWP7ej8w54P3f4Vh+JJ6mg9i7x+vePIkwQCzKglucOonrY1+kXM3CPg
W6uWQyeqtzyJlcGkBwRofFrUOIq8N+iO4sEUAqrU7UVzUb1qH0jrj6Cg+vpu
0Xge1XNb4vbdu9NQf15dSShD9W72uTiU8ZiAWJ1Qi7+dOOb8FNspkUnA3gWt
ntAOHPvCaRjSnEXAVGlGaCXqH387JuNiDgHRF4Tl3qP+ihZ9f/31U/T8lbvG
+1D/Fcdwd54uJiBPfYbDC/Wr2tQD6YUXBCSsVz3h2Yxj9ToS19JLCTjCrRfi
jPqb8Wef5MdyAujN20s10TyY9zzrfasa9ZOPuNTjWhy70T7esqsW3esc2psc
anCMa6f7lp63BOS3qplIVaP+mAxq2taA3of3i915jWPZGmtEmxoJYFSftVCu
xLHdj+I93JtRvh8oBDLRPKr4ubmBrxXNN8b1G57lOKZhliVcSSWgfNnLhRPN
r/8BUBxizw==
       "]]}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->True,
   PlotRange->Automatic,
   PlotRangeClipping->True], TraditionalForm]], "Output",
 CellChangeTimes->{3.529855415972924*^9}],

Cell[BoxData[
 FormBox[
  GraphicsBox[{
    {Hue[0.67, 0.6, 0.6], PointBox[CompressedData["
1:eJwVkHs4lIkCh5kMNk/ssEx2U5vl5KlWtNlOdtuf1c2uXBpKF+uwlayyybS5
RS2ONYwt5XLcTlJsNsMylMKoycwsk2Z831w/tzHR4YkZXRW2TvvH+7x/vs/z
rvz+OOswzcTE5NA7/va/uhMlZ1fFgNo5lSQT3sMOsrSoyikDv70fa/ZknQiJ
rbtfR3r+ireXPpuqzxODkOY0nHcuAu9cdgTZJcHjjmhd/NpyKGKC2etdemAa
5sxxK6/CQUFYml7ci/wgRvXSqiuYsU39R4pWig/9TNota2qxT1rSd8SyDzU+
RvJVXR3UXT8VPeY+gOem4amJhnqcfY+WnTojg8Czj67lN+Kwa0i+C18OF/uI
c4fCm/Bq5sXzH3L6EVugLB8cboZtKnPHjD+BpiUB10KjWnCVeXGQQScxy+lu
uf+wFQsB9ORaksSX5l/e3Rp9E5b352TMEgUyMvh9HRNtMFgddbkTrETPm9XU
hqO34SGO3M1YrIJ1avWj+ul2uO/WRrh2qBA66/jMJb4Te8xZeRERapSxC95W
PBXgwCneS9m8GjOJZtm53nfwwIF3Y/kvGth97sMVz9+BsfJXY8WYBp8/T71A
67wLjtxWejFUi33NN//zVboQzwqWRkk6tTgd/+y/KbiHzoE34a72FC65r6u5
YdKNwCL2E/0OCsKp2OtP73ZD8GFSReEpCuO/1za5Z4rgtPOvhsAKChY/6G/G
bhVjrnFs2SetFFavWi6opUtAtXm+zO2iEDi+r1svlkD6T2PzZQGFE1eKepfn
/Im+Jh+Plg4KhVH98v3f9OBG6nSetJnCzRVL1MWLe5E3Z7fzWCUF7ZDfECHt
BSvdeHvsXX+hPOuhdb4UNsxFQZF7KKzYf2fy28D7WHK8QtHoRcF36YIx26YP
P6f25NXZUIhWbXwplPdhfvZo8rb/acEpZC+8KXiAYz9yRKpWLepZjbQvQmSo
/ZjjbvxZiyM5P2V9sViOlDWfFh3arkXykYu/aGPlaAiKX8ZgaJG7vSk3USqH
84/9l2lyDSpcZfn2a/uxflFOadC7/zyz6fN8bj+mmD0mXsEadD1cXLhruh+6
S+zyFHMN5EK3EmMAAf20gckVqjF6eXtZfgOBFkZZ9uEENZ6dPVS5xoZE8O8Z
Ntlr1aBHZlT1HCchiVsXG6RTwQFVV47ISZhxMyJqilRw+3u3pwJ7o89wt+1S
YdNfA3VXChQIbz0gGH6rhP/g6/qvnypQ6T9QmtahRHg7848RlhInxcyQnHQl
4sq8+Gl8JQZ5dg4xG5RITw658dEHKsTtFnJfEQqc33vi1q2TKiynnfnah6tA
9cZzHWFKFSItDHa+3yrAd+B1vfBSI2XKMabaUQHRi17hxWI1NhG919e8JKFS
TIg8Z9UIvTBuPSEjMcE375GFaSCR7zSalpKYu+ByP65Ng7DwxKshkSSsEnxl
Vo5a+PrxN8+tJeG0K5KoS9ai5vpel/RFJNw90pU7KC3uRUx6m+oI+NhUaMa9
Kdzl1OcESggEG24NZJZT8PutiyhrJPB9n3p45QIFsuBsLC+fALv+xWhX+ADG
tqjsJ5MI/DvPbvy7zgF4VAsT+McIlMR6Tsw7DSJ4mZnnllMErn0T9Lg0fRBT
HmZlHDaB225xho0jg6AHWqacSiAgtch7osQQ9KF2uzZmERh6dO05u2oIHsrj
SXu4BAwi8SzDdBiL9n06SOMRMKkZm2uMGsZ89mjGo3YCtlm0NwHCYdgWSz4r
kRP45ODHplPOI6AHsAMMkwQ2+H5llps5AlrmieI5CxLbVoZbuI2NgFXYUG7v
QWKNyWb6QWcd4smrVZ1BJFqMUvNmlg6no8fk5rkkNo/stzTN1MFtMjl9fQcJ
8YPJ94L5OqzoEAfumSERJEiyuvRQB5EhzdTbRQE1z2KJwW4UDip5XVGCApGV
xdabt46ifbaA1d6mwCTX9X3uyVEcvLf+tcBKiYTTLYyBq6M4/SRrb1KMEnNH
t9itVo6iaYzSj3QpkXmA+CCZrgfHIfl21DIVrPyjHP700sO6bVrkV6ZCofcM
kxmtR0XasL+IoYbT6jOO0cV6yN9yHVdlqVHraP1Rq1gPEStpttRZg/8D3HI3
ZQ==
      "]]}, {{}, {}, 
     {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwd13k8VF8fB3CpX0VKlDWUJETSQlr4INFmifbdrkjSoqKNEmVLEiVbkiSV
Em32ZcaaGWNm7h1ikiIhVEL1nHn+8nq/3HvPMud8F3XnIw5u4mJiYpnjxMRE
f88o5TxQzV9vukio3+fhWo4TTRyZojgt6N+PF5/kVQ6ly7fvXYtbDq1JB50f
Hy2HV0tgz9E4a3REaQ5Kny9HrETvdLe47bjxaHqt7u1yZCbPv7cvzgMnVPt7
XzWWI7fbwnt3nD8C+XnrpKwqsHx64/udcVfQuUDtuodFJfzsjvLSxq6A1XFM
1nVjJbIjZdq+OodCZZJDqseWSsyd6tB3dlEYnGUu1wa7V2KaJHtqGuMqqI7r
t6WvVaJjPGdD90gEWtbQs49wKsEI8HZ+tygG3GeTelSPVOHEEV+BTWw8nHS2
Xit4wYD/Qe07ayzSkDtHRurIxho4NFpGuqtkYsaaO59rJtXDa//U3uATObi+
XO9ufdR7VMs5FcsvzkVeVkeq7X8s/ImUKf7ingdWifoIHcOGmQZP+C6iACOH
ZVvWTuLgavPM2TvPvYHaOWrYc0MzBhMfXD+aUwhqRVtXTQgXS8x507SbihEU
uvLl/VAelp43Yv3hF6Mq41r9hnAelr2Lvcn+UAztqc6nv0XxYLRis8r5r8XY
fflWsO4tHlYurtZpFi/BugfLt169z4OF+lvLS4tL8Mvk83pWCQ8O41LOtEWW
wJ+vmRQyzINfiUdn/PpSlMarquzbx0fCoyXp3+xKUWumvTjdiY/im3+cLLaV
4lDey9WfXfmQPnSj5atzKWL1pQ1dvfjIli1qMg0sxV4vwe0Vp/jodJEv68gp
hT23smTzdT52TaxINphZBqH3LsllJXxc7I/ae1m5DFZ/3PoVyvnIpHbNouaU
YXj6nXfDlXz8yumPC15YBnmpCVLPavmI3a4S2WxVhvbJb33GuHzUZx4LDDxd
huJL7mJ/v/FhvnHuTkZrGcLdJGyzFChkqM+597SjDC2jOiX/lChIDat+i+8u
Q2XppjMOKhSa05UuHvxZhmz7yYV9cyh4/53+UHJaOW4qZ1yRWEAh7tm/4U0m
5Byfb5Z4vYpCt0JLfOOdcni2vVs5bS8Fu15K+Cq1HO26x02l9lN4Uc7TS3tQ
jmafz5cnOVE4f7Sp2O95OQSWkm0/XCnMrKn5MqOmHPMPPSh55k0B514bbxsp
BytuiXlGAIXYjls8akcFvFSK04NuUZjjIa5qsb8CXNWcG5sSKDzq8nZ66FaB
wo8xtjPuUCjtNe/296uAu5mYzq0kCv3DX0dnhldg2uDoiVP3KdhMNZttV1QB
T8fR1NRcCryoLJeXFRVQuNJ10PgFBVdZuUzV2gpoNwSMq82jEKjQbdDDq8DJ
fdkpXQUUstRjLcIGKiD21ufB30IKE42+uJVpViJOiacSXU0hpsAha4FeJVIX
1o8fqqGguupd7/UllRiSiynaUkdhmVmM/wGQ5wVi+VLvKThvWB32Z0clQkxX
PN3GoVC4Lzp7eXgl3G2K5zt+oLChfaQ/KaYSx8Zw9VIbBY6Lm+HEhEoMbNVh
5rZT6PFcWdh0vxInomXo/zooKB/raDhaVIn4+X57Qr5Q8L9iPPhogMznC2vN
934KxkZdJj+GK3EX8kNfv1MY6UgINf1XCY/Q1PSPAxTOWYyqNk6pgoZT2uTq
IQqX/rxb93NeFTh5M7/6DVOwyva5Ad0qpFtty976m8Lk3bNbQxdXwWKDpp/h
CIVrBReOzTKtgoRaxtSeUbIfx9YkmW2vwjVXmVjDfxS2zB36Era3CkctYsPE
xGjINaYvZbtU4XfEilAGcbz+RKabbxVajW3u2YnTSO5iDF4Nq8LNM2FOK/+j
kb3fdj3nTRUu2L8PC5Kk4TPt3w21UjL+Zd/dalNoLHr3pNWDUYXXalZm+cS5
yjLHR5qq0NMFqzYpGq847KTZvVXYcej4NmlpGmcuXeryHKrCE8nCphvEq5ca
LssdqUKxdYGX3HQaxdFxTMtJDKivjemUlqFRuXHn0ME5DMh5d2oIZWmEjkjg
+XwG/v14GWA3g8aGh6/DxvQYuGj15/sr4rqJKrOjVjCQcPzj6pCZNCLzag/y
wIC3ufiML8T2rmdfqFsx0PosVtJajkZTSev6Fw4MVOerOY0Rx/lGxf7ZwUCO
YVKJgzyN7bPNPljtZyBifqpVBjEVmHqc78VA/mGbEksFGu3LXZP/XmKgscw8
ar8ijTaW9Un6KgMGs7b+ui3yYV2bgmgGbLwzvrBFvjcwcjSRzLcw/aCpErFp
c6NtGgP3HqePHCH+wH+VqZvJwLWjoQrJxK3SF7d9es7Az01Zv36InOW6sPQV
A2HN5m9VlYnXrpuQXMTAkIKv9BqR23TpgAoGNjR8/O5G3BIgnbujhoHzoZtc
Q4gF8oOhho0MCBJ3eaSL/Kx5vyyXgYAlRaNFIm96bdQnYODNKr25fJE/351a
K2Sg8OQWuk/koIsdmV8YoL2/ak2YRazq9uZyLwMqfvcny4tcsC7GeYiBp+nL
z2kS0456BzHCgHPP9rAlxFSvtJmKGBPKI28XmIgcNij/+z8mLK3VXNaKPI/7
jTOFiW1yy5ZuFLnodXmuDBM4nH/HVuRdSXeiFJj4oOZ8107kHxf9vFWZOFvc
a/z//0e7rV+vQZ7PUTyygZivu37OfB0mkiUSzSyJeZV6v8QXMdFWYJizSmSn
6fUfljFxxf1egYHIY4Ppb1cyUXvh3gENkW9xAxLMmJi54kfGDJGXvHE4acWE
a/CO8HEi1yXpOG5igisXKfON7A/PM0jMwIEJGc/ty5uJuePduVI7mCjIPfn3
LXFz0vqcrr1MdPbmeKSKvGLh5UoXJhI5ZceCRT4ytOTCESZ+WsT5m4ksyZPY
e4Lsz89W31ki33/TtiKAie3ZytOHyHng0EGRgyFMDCSFrUkibjrp7vY+nInY
7zZ80Xlqktmw+nEME/J+WooQ2Vqm2z2JiRRDy9hmch6bhEPFa9KZuO+uxk4S
+Szv1pwsJjxSVhW6ErOfJ6+l8pi49fd4Yhc57yzbYJX8N2Q+3IXXM0Xuch+8
UcLE9Dkuhm4iz9ZPtalj4qnVpSguuS+s1zKnFrCZ2Bohax1B3Lj1h+0kPhl/
pvMzM+L3196OFXcwEXXmemoKuX8NvzbsWDZG5tvnO2ktua9Dpnd0osWrUT6b
EdBG7rdyyNeRr5OrkS8ptDhN7D4z/O49uWqSf5/WpZH48NegTii7qBpZs/SH
OCSezDul+sLHsBpngiVV9xBvKDp8uXpVNeaIe6R+IPEnzmaadtC6aixJn17S
Oo2G/iG7w9+dqnH4gnfAOxK/tjxLNrXxrMYjzbBkfeIzw33SD32qccVVTfsu
iXeVIdG5BwKqwVE38/Ij8XHfvcZfDbHV8Ej35f+ZROKPwDEop6oaE9O6diSQ
+PpCI91Rsr4aQRt/2baPI+f10NA896Zq0KuU3s4nnv87tkq1vRpyI/vEskj8
LpRrlooYqUam7CWHmDEKvbY7bnktrIG4e6De1x8kn8dlelYtrQFH+e68KcQr
W4ZXaKysgeDj2gfaJJ+EeCUIKKsaHG7X3ryH5Bu1UEp9w4EamL4w2JXRS/Jz
yZ5s7Rs1OKeUu7uyk0Lu0AcJxYQayCbJaTd+oiCv7eIxKbkGfxPvt/NJvmuN
PDi3M6sGpaln53wSUvDZczL+XmkNbg5GdgtbKUT+irqkNlADve+SSYHNFOoX
lu2ZubkWkdXfpAzKKSxxsno9fnst+jgR10dKSX0Ty1QY3FML5aQBzbISCvtG
G1iNnrWQdDwfZFtE1ssQrIu6UIulFwd0bF9TmOb6c5nk01o801uqP/CE1D+3
daaKSdeBbR06nkPqE+q0kpuPXB3G9dYUWcaT+mKnxDvBrDpc+HHwQG4chVOK
Xd6vtOuQX5655coNCmlxmbV+FnV47P57iWIEhaEYrfBPJ+qwQL1bnnueQvw1
zSk1gjpYndFTaHKh0BaoLhGXVY9we6O96vMppIRu9+h8Wo9B5e6S2fMoOMVG
VBjl1yNpXe1ZlbkUhI9+X+SW1SOgVmucjBqFj1TjiGJLPQymphzolKPQaRzU
c0e6AZqtWc1r/iP1xZCwIfV4Ay6Lc2yahHwEGjnUbJv9Hks/ZiSmJvCxweio
YhP9Hk37x5SGx/PR9lrzuUFmI7Z8jkk+4ciDcdXypr/uLFQEpo6rzOBiyy7p
wqlz2Tj1yiB68mgzjP0Cf4z2sPF4uTJHe00zWmWsn+/OasLtWR186yQONho9
KLztz4GmiUnX+l9NCGyvjVE2J++5nWR8d27CkZn236b8aQaiYrOXd7Fh69+/
aJUYFwl7Gdq2n9nQo6L9Do3nwv7U3WG3T2x0Jzf+YkhwQXtp3kxoZ8Nd11E8
RJ4Lv2zHhukUG/sstimIGXDxIfFVLmrY2OS7x3zImYuS/ND9aTlsLKj1vClg
cJHruWL4+mk2tO/PTeLUcKFdRjvPOsXG/HOCjPp6LlZc9Z+ScZKNuQb2BcVN
XLibRikUHWNDOXY5nd7GRXscQ13Mhw2JPRPVDw9zodTx3KDQmY0vX+9l/9Hi
IfvxsUe6Nmx0lu/N+7mAB98LEmYtG9nouKtQ2LeQB/t+WaXoDWy02V1raF/K
wzYx7YTf1mzwnvuR6pSHKSvrOLQFm/Rz5saR23n4V/jAucWYjQdTPpSpXeFB
aqCvgKPJRvHIm4bJV3kYyO/tzJvHBtUVTw+QPmoyJ9v+lgYbUgyHwcoYHm4e
LOTuVWfD91Ll3CNJPMy4Hvjytwobxn8eXyzKI+/XNFV6zyTj9QbiQAcPlGbT
bbMJbLS37Ny44TMP0XWyTuvGszFSa7R9WTcPP+5Wu24WJ7/Po36fyf080DKv
53mIsRHt4Zr0dJSH8D6JK6ljLOxo2/hHbAYfn6s2cTf/YOFLo/KbVHM+fK9m
3U75zMJMBX3LGks+TH/n7i/sZMFsj3ndkDUfS8/Mu93yiYX4To8P1rZ8TC+W
ODKng4X1oy/G9+7iY66L9q38Nhay5tltWnmcj0Ru43tdioXmg84cF38+Cr22
dhzgsyD+5MS+iDN8NNlKe9/isbBrRaJv2wU+4h+XcCdzWZhi2xUbEsHHNlkT
1T9sFrz9g1tYGXxMnGZhrVtPxn8b5z72kI/9zZOu+9exUD4uq0/zMR+qkZo2
FbUsqIS/H3f6OR/Qoma51bBQl6I6f3YxH1wJVemXDBb0q/N9vPh85HVdX3C3
jIV+lZ5/46Uo8O/y/517zcIqnqLdr2kU9AstF4oTh9xYm9QtQ/quCM34kFfk
+5LJqxtJ3xeVGaN5o4Cs/9fm00kaFBiFEspvXrKQ1pg/YEz6OHn9fbB6zsK3
8A5zPVMK60I+L2zJZcF4ncz12eYUDmdQ204QNxQeWjTRmsI4M93tD5+xMPZI
1ZvtQOJw3K0pc56ysDUk+NPhQxTs3yyese4xCynmT5c5HaYgF/gkciCbhZ4x
QfAWXwpsCQXzJOKgY0ZzV52kEDHpxfrfj1jIOdC1b3IQBc3DFz6+y2Lh9yz5
nNFLpA+8liXtR2zJtfjTe4XCDbkmXy1i2ibxDofEye1/tBpvPmRh8io7XhqJ
u4N/a8UvZrLg+DNAKy6RgrtrX4MpcdKzzJNhyRTCDKa8//OAhWXa4nK+pC/U
b9PwPE/sJJe32eQZBZMFx/ujM1jIft+eYkD6wneb/tnsIB6+Nq1fI5+C9gNn
wRziSHHPSMl3FG6eTLiad5+F133KNVwGhQ9uEzS/pbOgOO7cLcla0keyxJRe
E5+UbXcxaSDrs7jpG0q81PDBWBrp+7rTpNfrEMdYSTI5PApln42Zo/fI7739
8M3JAgrt3AuPGoizzyzV9yF5LnH5t9gAYsnwuJEUkhd3XLS/40jsefd3JZv0
fbt07cUXEmsWF+1f0Ufh28rE/o9pLFxqnKvnTfKuUErDv4RYKLw8nETydJri
mz0pxEn/bbo+geT1ZY8j1jgTj8k/2buc5P1ks79r1hLv0pZdcIjUDQyNnlQd
YoWNvNKGyTR+/bh//GcqWd+eVVHipC5hJdz70krcdDhptyGpa5p/XSxgEEdH
uw7dJnWVnKvpoRTi3tSq4jpS151qWb8xktjm+YIIMRUaKl+iQ84SZ5dH7Fw6
m/SFnzJn+RBLNvdrus+lYfx0vuQBYs/PjgPxmqQvupO02ZG4cvhlYY02Deax
9O/WxPMkla/91SXPX6jpMiEOnnV2++JFNC5fzFhhSCzUa9NwXUL6xudVbQuJ
zUzX9McZkrqoprhVizjJLuMt05hG79jcpRrEYwckwsZW0Tg/dK5lNvEuP++t
i0BjocC1RZX4VXCDurMFDetdrktEVry5pDd2LQ2JGCmByCczbr6uWkfqQo0e
vuj9pvzhkJGNNFz0r+uKvr+EudtxoR2NcJuw96Lxo6nC2QccaFglOjFE8+v7
qt4Ts5WGyYsMWdH8bf9cKqjYQWNvTkeOaH3Z075cGt5Nw0jjQaJo/ZJzNm7W
3U++r3uadhDt1+Ic1X3ONHr2fzq0X7RfFjLd0W40sqgTGw6L9mvL8ZdlnjQe
G704HUgc5MYN+ulFIylaZzicuP3kSjudIzS0ruyrSBLtV+jdWXv8aGxdRbc9
E+1XgtiXyBNkffZaGypF+/O28sJQAI0DxesnDBIX1OnYaJ2nYTqot2UKOT/y
H8KVdgWR8Vf+7ptHzB7nmFsUSkPay0liD/GSGS/PDVyjYeE4MeQ0cfQ8pY2a
UTQs7TI2xRNvsv7w8epNGsssO/L4xMz63SevxNP4L8VEfITYahtv8qU7pK4N
TvWeRe6DuVvjwrOpNKqP9SQ4ERsFl/kfeUwjeOqyjaL79kLSTNL7KQ3bNWGB
88h9XBzzNtHzOQ3dwQyGLbFeal6J0ysaCjNin2QQqxc9mLKlgsZt4/BsF3L/
k6zmJdkzaIzFnsu8SaxSn2JgU0PDWaX9DYNYviVhq1UjDU0JS60lJH5IjV5L
Nm6h4SMVPkGWxJ+fRr5LVX/Q0E51PvUfiWfHCnsqlIZpTH6Yn7+OuH/twR3y
ozR0Kq9NiiD+usX5rPQ4AQyiu5oVSTxs89tSNW6qAEcLB7VMSDytyTHe3TlP
gPa+sYUFJB6vM8r7JtQSwGu/IaRyyO//bvGFDwsEuJn3cKcTcXHtgvs8AwFy
gg5lTnvCwstulb7q1QLo1axMOUHie8r8ccFPtghw3iUy6jjJF8blKYp3twuw
3+TtOC5xo5NZztVdArR9i/Na+YLkz7vn+W4HBBiRvThvQh4LLjP+LlLxJt+/
vVftHsk/88VGBaHBArjqSNRMI/krmz9k5JwrwLdx0SUJJSTe+8fW2uUJcNL9
h4VyKQstM5c5mxQIoLPuX/Yd4mn2xyIUCwXgekmvSCX58Ujl94/1TAG02wI1
8ivI+XneG72qXYBt1jnsqdXk/IV/6Z45vQWdJopaLiRfu14o/GM+owWBSmvk
pzaxIH08dvoR+RaMdm2dUEDsvtvMiKnSgl6nOGp6M6k3FsRfPKvTAscw7b31
pB7wrbJW7LRogXl3x/xTrSxoiz9Ym3+iBT726nODukn+93dL2Um3YHWMMHn9
f2w4t285bd/aggnOhm/UJrKxcKOlg3V7C07XxkUOEZeqaUww/NwCSwvj3WmT
2eipaPOYPtgC7GYFTJRiw2zmnkVVkq14ymnf0S9L6sWnDoXLVraCFenga0Xq
sxXdEEgntOLXegl90i0jrf2J2333DxhZZqeSFEaefzJAmyi2IQA+fdF8Uj8u
Wc1wOtkG78axzfbaTVh1yn73CUEbrmw09Yk43wT/izeelaxqR/Tv8k9u75sg
Z7vCWz+hHUenHVaXU+VAoeVh3dW/7XgSNZXl4cXByKL1A7P3CvFskcDRNoeD
4HCWVeU+IdSrr7gdfMKBVNfuO94HhLi+POdT8FMO1NJ8LF+5CLEur31Rfi4H
FjNj4xwOCTEvzKhdPp+Dq8OtKy/7C8FuSXR5XcTBrJLjwd3XhWhnrNkX9Z6D
dNW/3OgbQuzteRZ6q5ED/TNX9JbfFGLw4mz5ZBYHZktvc4LjhYh+mbEju4kD
1/tFWqrJ5P1bJseKeBxkh0nW2mULsfBA4cuqNg4MO2/M+flYCOcOt+6Kdg4K
LVRPJD4Rwt/6Z0CZkIP3o4vUunKFmGB3J/9tBwdDh7f6Br0SImF/TOiDLxyc
rW4t134jRESV0r2ULg4mankqNbwV4laetlJCNweKbWdKZxWT9VW+kA/t4SB1
9QSF0hIh7jQm3T33jYMFCRFenmVCyBh2nzvey8HzH/LF0yqEaA20fX2wj4PV
Dikz8yqFsLwfbLevn4P/AXUuhC0=
       "]]}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->True,
   PlotRange->Automatic,
   PlotRangeClipping->True], TraditionalForm]], "Output",
 CellChangeTimes->{3.5298554159859037`*^9}],

Cell[BoxData[
 FormBox[
  GraphicsBox[{
    {Hue[0.67, 0.6, 0.6], PointBox[CompressedData["
1:eJwVkHk81AkDh7GOYrEl1/rQkkplHe9Wovi2IqnEat9ssoVaTVNtoVZWl2Mj
RxkrbehtIt6UcoxWzobc4xjHGMPcjGYymF/SJNK2fzyf54/nv8cy9Iz/L2oq
KirHPvOvjzRFtV5dS0JMSdalbVoD8Oq/c4tqHoeUtpDmdfQBRD3771yw4w1Q
tlzjeZJY6GMkPU23uoUe6kFXx08sTNSGCc/a5qDq/tV+WtEgVAOsrtvkUOGm
wzjZe4qNNN9leSbUfNgEfuH5rfUQvt6lUrOkoBD24bEkc9kQCrZP978vKsKR
R9Mk3QEOHJ35cunTYtC1vKpYJcOod+zS4NBK4J6mzdbLHYG14eGbx4LK8KB0
W0JwKRdkCiuHyy9H3on/n/uljocyXZ+HP4ZUgGp2XmrRxYfyelNF5+gzRCh/
rw6cF2Cb5rYGj7BKLHqblTXfFiIujtZVK30OZVbH11EOIrQvrh/eeLIatjzT
WbtuEfRi8saLJ2vgX5dsIjwrxo9K0xnrs3XYa9waXb58FNmRlE+5b+qxOL5m
X3fDKBRR6teSXejQtmjj3b40BoPN21Nb5uk4EZHECXKWYPPbmAy1ugbMVPxn
g9GcBAfLK/9yu9yIDF+n0CMd47h4duZ/v+MlKN/vqCpPe4V7dvYFf6s0wfBo
1ozRbika5eTHbxqa4HdXndVuLIPkUWGZXXwzauRfNJziyqB1QlxJ9miBu3xs
3rvkNdavtagv1GhF5QNubFXcBPZJDjaJW1pxaourio+bHOH5tzosktrwrN3D
U7wgR2ZILzPQux1WyfMpj7snUblSl52l3YFEC/lY6Z0pcHi7eH2MDnjFuhUy
d09jISdhVC+Ngb32CTGyuWmsDKTLdu/rhFj/J38tTwXcTRamr+l3YeUaE7+N
yQqEDTq9a2R2ockyRGLdqcD1zMiFRUo3SEl18ntLCRT7l6ht3d+DnXZuL7c4
EziedD5hqzYTb5OuXHALIxB9/M9EDpmJImVwvk86geSdZclRDCaascZ+3VMC
uat70gxte5F+rGr+YzuBJ+qT6bTUXgj5MXHhXAIvRrUzf5jshe2E+lbKawLM
Rpvb0z59ULL/2l+rICC6vzM77WkfvGPnFr+cIjBz9djdDfr9iPxYsMxRQkAj
OI7afqYfQSXgdXIIGIGaf5zZj9obO8dz6ARs/t3tOICaNvPed/cJOH8cKcqn
DMBvfuUNrzgCe7hzxd+/GcBXZFEw82cCQTXGpQJ/FmirsvbCicDp7E20SzQW
MjwSHxZoEbgcvf9vsxWDyPv2faJDhwLpP4VXVZ0bRIRXhuVjigJ5TjdrA1iD
uJipDJbsU4Bm9OTF7CY2HJacsrdUU6B5tqPxzyw27LxdGr1J0xgckDY7Ktk4
6e4ZQG6ZgpSm2d4TMATnS0P0OLMpfMiw7jz9fAiLyavGSn0moRPh3qNjygEj
acFOliSH+Q/BfUXRHAgi1qpvr56AncNlltfw5162tvg86zW26+cOSVyGQbrR
e3dkQga/qaqR+Jxh2Dt0+mppyxDaxeZbLgyjpql1dZy6FJHFs6IXQSMYKHX6
denyV/gjxUDyc90Ihkx+MxUZjuM22VE6b84FtbRbzepLCR56+07cucyFi2ao
vabuGKptTk85CbjY5SaYSDMcBUMrhWCBBwNSLOnRJxF44w/fRlJ5qB0gr9gx
K8RUc4tymSofKRTJme/YAqgUjH0oCeGDUcvv7MvmYXmC2qJPIx+HrhUceB4x
glVHv1GVWwmw1+q4l64/Bxvd3dST4wWwfZm78GAzG56WQVo2YwLUf6x+FebA
wgYVV42jVkKwbA7/EbCjHxXTDM1yfyFG772w5P/aC1dB4BLVeCEOGme+Ep3u
QUu3bKkfTYgmL3f27M1O+NZf0Lk3KgSHUvn+Tms72E+0dKcMRDgfsrSRHtWK
4LtZeq4eIjzWDxUGZDZDlrr6q9RzIhgdUD1xSPYSERcrlo08EMF9JixlHbkR
H07uMFjPEiHN751/dngD4g/1rYjWEGP3lbkLUjEdOntCjNo2iSFO/GbwfSYd
mS4KY+MwMUpzFg/0pdJhvv6KaViWGBcGY1hyGR2Fpnpmz1rEWHd0z+G23Ab8
A2P4Pwc=
      "]]}, {{}, {}, 
     {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwV13lcjF0bB3AqIZRImwqpyBZFhfhRJCWRpYQkS5tQKpV4pN1TIUqLp1Ja
aJEKZdp3rTP3LCWVZiIp1J2ENu95/5rP93POde7r7GeW2V4yPyc0bdo0z+nT
pv3/11suO03x9d7tyQZvz1GRXLhzuJKlUSug/k9sil0UF3IBscn/RumgLv+X
ptBDLpw6fb66RO1BIuuSzp5oLh7M/j7/XJQFhO2rv8+P42LfUIzM+hArnCh1
TmISC7fuVhq/dhK0xsuUyEdcbKaaI+tyTkNWe+EjjXgu0hPUkq2j7MArtd2f
/JgLkfs+alUSDtAKPhMQlsSFTSD1VD3EEfoPauZdT+ZCxvlGzsg1Z3zdfn3h
hRQuDBQHrKtzXHBm1sCvkKdc5PbrXzgedRVu72KMq3NJ/Mdva8UPeULPoD5T
JJ8L8a7owTIJL6zYNGy79yUXjqzvrmoh3ljy+XGD4DUXywpivehr19Gvffve
kWIu6vS9KitzfJFtnCK89S0XCaX7VRQlbsFD4WNOSz0XHnoqAR7OtyD7z4U8
h0YuVLRbDNVX+eGszK2CrBYubqqr1ocn+yPbKKLSnceFznwW81hUEGJvpH2Y
8YkLVzOXtqSJIDiHXDsm95mLzHDJ7gHbYMTETclofeFCeZ754HWNEDxZ1mHr
9Y3kL8ael1R3GwM/5dYdHOXiozDXuH8sDJpb/wasnc3DEgN3c63T4SgsY+5I
mcOD1a1FVj614TjxIGv7cnEeWqYddRB/cAcFpgIhzYU8dN4YMIiJugtBz670
WEUe6q5dsC3WiABnxG9ztBYP5+P0evm2EZh1ijXPTZsHYcZcR9GoCDhPlS07
upmH7WOZLmYTEbhz5eLMteAhz/PbTX7dfZQ2/ChXNeEh3v1iwgzbSKyess7S
P8uD+6XLHaYPojH/40TTijgehOjWcvGqaEi0smeGxfNw1xVpLT+icW9e6dNf
j3nIdBd3PXgoBppLNql0pvPw8VrmzCMLYqH7eWZf3yseDgf3aZ64GwfxDdLG
whweNB/bhDiGxuP6naUrdixoRZly3cXVRfFQrKj1PSjdiv1PNA5/HYhHfYHI
Jnv5VjikTSldNEnAbeufSSnKrWAe3d+uOyMRa/KOh9potWLWya+X9KoSUT5X
UZx3pBVXHVbGGegnoafBy5QR3wrp9yNqmceS0HuU99I6uRWv9pXnSrkkYTVC
1s5Mb8WIhlX9p4QkpKWqGdm9aIXr6L9/giaTgP8cpM9WteKi36BFY0EyNkkY
u5kNtOLsf6+kjmikYHbMoq7Pem0wZ+0KP6+Qjt1SdxKGvrfh9o9K4YZN6YjV
Ujh58GcbKhcZeGmYpYPi1zW/Hm/DRqudZ3/7pkP1zfG3kbPeYVHPti23P6dj
b5TIr0jld2gb0e7Nyn2KMxk1c1It3sFaVn37T6MMmMtK3hivfgenU/O++7ln
Y0zCQ009rx31i06XSW/IRdfnJ027UjswGS5Z1nf+JZRKKn8NSnzAjuVtguKw
AsRe3xz4+F03bvOklhy7wUCJseRSwToBfjxKu+eSXQKXHZduKxX0QHNnm/hK
Thk2ybsa3dL5BNdyu97ovRWoHV29eVlhL3aaKB+r66qEy+8gXzGbPjz4+LCt
3bIawQ7V+zL7vuBqkO6PjOEaXLAd2CUSMgC+ztmEKf86OH9/NG+z5De0/DK2
3DjxFofrTwmPOXyHafmJzJX3G/DxiWyNTcUgzGLV502TaMLc1WsOabgMod1L
7tzFRU345djHv3plCGePzS7uWNyEy5sdJErdh+Ap++VC4comXO9PM9jrPYSk
qPRGV/0mmLW8Ht/nN4SRiBWhn9ybkCRYOU0ycgjR/6rOaehoQq3qRoXA10Po
9lk2O+pZMzKnu3PGfw0hMdjCrjenGWP9Kyw7/wzh9IOwau3Xzdi3vsCqeHwI
gow/vq2VzXgW3v3O8+8QetpZY7Kdzfg97JrdKUqjV/fW1ziJFixq7Va2XkTj
64ig5bFbC6zua3Zt2kAjc5rcumHvFqRVRYj2atJwnmsWqu/bAmXx6bEPNtL4
vpxh1BPWAunFaX+/6tAYNL9frpLWgtzfwS9vb6cxnGOQl/auBVoGdkZeJjR8
tM0bji5hQoHSnLHkLI1ruYFlRspMBNaWjvqeo+G9jvFyiyoTuoypqz3naXip
qSQqrWbilU7TiiQHGh7So2692kwoi1quF79E4/JojJLHfiZ2Pg3wS/SiceYV
3yXyOhMP/nuzKSOMhq2mtF3QTSZMvZv3fAgnzjY+4eXHxENRoQ7JuzRsUvMM
T4YwIb91dYhrBI0TUQGLVSOZEI8zHVj2kMZRD/Xq/EwmqtYv9dBNpGGs7SLL
ec/Eeqfn9M4XpP6sjtszPjBx9G2O2fFcGpfaDSe1BUz0bTnc5ZpH48ENBUFM
HxP2R57+++gljc6a2oxTP5mIulci115I8rdQwoA4C7bvX8RPltO4pR7yQmEB
C3vMrjP/VJD48R/L9y9iYbdL+eqRShqFCW9nvVjMwnmhsRM91TSE+twoj5Us
4K9mfMZbGpGeDeeE9FkYvxO2rYdJI91Yu01rNwsHTilxylk03ig8Nj5nxIJP
sWdYPEWjq8xDo24/C/mHB8+bc2isnK38J+w4C6ERlU+TW0n9GM9QOXcWOh8+
lvnaSaPRqWfK2JOFLVHPMjK7aHzYtt/F5xox2/iE0wcaIoLlRz/4sjCnnCvZ
003DdFXLkpRwFoIsqUpGD43uN6p569NZYDLpYckvNBYIPxtYncFCS1eZy1Ni
A5N1KiuyWejfqm2Ifhqp77UjFfNZ6FP58/3cAA2niT1XxcpYMBDfopvwjcbo
NsetH1tZsHpV6lxOk/wDv1/50M7C1xnxatuHaRxrds1s72RB3v/j+kLiIutr
iqweFo4opes++0HD95/Qv8WDLGxcv0Dc8yeNueXZlQ9nUjC/MRVT8ZvG9tla
ExFiFJpifpqo/CHze/D1xvB5FGbke5v4E1P80hT/hRQc1Qw1MUbjoRAryGUJ
BVo0LjpxnEad8ZGyC8oUZC0KckaIxyLe/bZTpfDW55HUngmyfpYLHKxXU9i+
q7Okl1h51w8TEx0K96Yfeig5ReNw6FV/wy0U1ru+mmlJHMAZL9q5jQKktCoe
EfedFVmna0DB+3W61LK/NOSygs9rGZL6S7Zm2hAb/5ybsG4vhagAubB44qyA
RfNVzSjUBdPqi6YNo6spxmipOYW8qqBPpsTzpZV8Fx+h8Gpb8hd/YtdUtWHJ
4xREotLqvxHr1upwps5TcDHL7H82fRj+qh8kBhwozM63luQSs/wCTVovUOAs
Nrg7Sey4g1vx3JVC6KxfkXuFhvEq3mcqzp3Cs5muCs7EwpPLtwR7UgiYch8L
J35U6PrC5gYF7YErRU3EfTLy3/b5UhC3VkjsJ97kUb5ysz+F2opn7TOEh3GL
Y39WNYjCvr2/nJYQN2vOT5S8TeJvtR/SIZa/9/r9ZCgFgfXaCFPi84PWMv13
KHhteKtmS5xrKnqIF0FBLCtI0p34b0ZWeEUkGa9P+uaBxPvEjtRnR1NQPMDt
iySOtp+YERdHxqdCjptM/LEmeWdQPIVHbgPSOcTrVU2uX3lMQW9KNfsNsY/f
cMGpJxSMjqXEVhK/5ceMmKRRqOre3VFPLL1j53rdZxSUDn50ZRLbxvc5qWRR
2FB30ppDnD1xJ21+DgXb6jvxPOJxK52eiVwKqeqWWq3ERoVdSl9eUojRiVT4
f/kDmUArbgGFwPuqVmzibve1UeUMCtcN+gaaiddyOKysEgo2z4tYdcTemj7z
YsspDCfemlNOXHN3+d7AKgq9BvIxr4kXDNb7u9aSctWz3pnEp0xdy6zryfxq
73yeQJyZITdh3EThcl2Ezj3i37PLdXSYFKjwLQt9iXfZ219Zzqagnyivf4n4
Xo3Ecwke2V/ui6qOE6/ys1bre0/idT+XaRB78GfYcrooxAoHbJUhrkTWf2V8
CiX7GmZNkvk/MTEuFfOZwiKznthK4nSr5AMB/RQm/tw49YR4pMA41OUbhVOi
Fh5+xGHuMcLGPygcO3c4Vo+4nb0D2qMUJksN06SJ1TT7vJX/kHxCJUQHyXot
+a49PDZFoau+oiCWWMy0a+3n6WxsPHl97kVii4wAB7YIGxGPs1+AeMiO050h
xkZapdxQJ9kfejXXFKLnsaH3Xi3kGXGIynJL//lsPJM9eMWdWJnv0nJCmg3Z
krQtosSHrSRK5imzERKkdH8h2Z/54t88NquycdP8tWIT2d8LK+s1zq1kY6zb
LdOfmFod+LhoHRvlsdN53yZpHJic8HfaysbwzF9H08n5kfOifdvD7WzMyd4z
zZxY4nzBaMVONvIbJbP+kPOnufmKvbwRG12ytfMNiE0T+03eHib5caumisl5
lnm4TuSnBRtS1l/LLYjnzk4tXnqcjS+1IsGD5DxsdDmt4Xmajc9GzQpyxMYG
bQtWXGTjYPk+u8OjNAx7q9/5B7Hh1DDTw4mct6mxyRE5t9mY/89920/kfJ5h
5mvSEcZGgJ/dwZPEVa/1irUesGFxg7VpzxAN/ZC8REEiG/O+rpIT/k4Dax7b
7XjDhuHu7D2Lyf0Q331j6YViNtRaNCzc+2hMPTjx7mEZG6KWnfZNn2mUTMqY
DNawURtRF3G1l4ZeS/i6eDYbNqWHpueT+0fX1efn+Fc25opn7W8i99mrR1fc
6EE2pEWn10sRb6p1/NE7zMaCkcIQqw4amgrHaNZvNm7J+/K622msrdn0LU2E
g7l0cGobuR+V5Qc/HlLg4Giu0Gp3cr/OqTjNyTDh4OobtoYtua9vf7U89Hg/
B2t3mpu7kvt8lswBKuogB7GG2oyb5L4XvbCd6WvBgbb8vzbRpTSmSS9utDjD
QUHCjIlcBo0RB06l8DUOglLDQyPJe6JLck/e8WcccIWTOuYnkfdDhNvzu1kc
vFgyvIlF3idCC5MyqnM4+E+mO/ROAg1VqclkjdccTO05sFH0PxoXpPMeCFdx
MCBWfPE9ed/8kV/ikdHJQb9QWLEqeS9JqfzSHZ/PxZNT9djqQebvicqm9VJc
xLruzmW4kfFTNd9wTob8DzY6oL7lCo3jalnqLYpcCPnKq224TCNp5Rn55FVc
vIytOT3Tkcb6tS3jxru4OOBxT/HISRom2mklsVe52GIjkXpVn8bvxl/FZd5c
uM0btS7fQePJGaPi3utcbHd8pTMbZH7vfmFo+nMhGzPd+P5WGrkDqwvr73Jh
nxo2dpe8P2Uf5+SNkf/pu2L/Y2uq0fg45026VQdp7z2jYsEs8t7kN0bI7+RB
7PyObZY1QzD6rm0ptouH/BfNy8yrhiA1nqg4ZshDBG56GVcMIUPKPf3dPh56
U4r8tEvIe36PYkm0JQ9sy+J5Ey+HsDnb+Yv0ZR4u+sN4W8oQRr3FdyxM4MH0
3vDcuFtDuCR14NucSR4S18hsb9cdwqpG+8iOulY8Kd7bFS42iLQ5HyqVgtrQ
/qz/dVL4N/Sx5BmPd74DM1b+kdbPAQwpfP0rPLcdJaVCLvmL+vFmUL6hta4d
TeXd+wu39mHfng89tyPfY+3Yp9galV4kqk33e364A5nMG+orVD+iILSvX2p+
J8IXv++qXS1A9NVzicfed6KrJndIcm43NvejQyKmC57jmxoUt3Qgif/8XMr5
DwjNdbzWpt2GvufD77fJdsMz3eX080AOujX16k57dKP43v14xVMsbPU8cNy9
oxsKKxpnvUxthI3/09Dhzm4cFjMZjrrXCP+7QiWXPnSj70TxlQCfRjSn5y91
FHRD99cci2DzRti+k+217uvGq6gsY7W/Dbi9hX/ZaKQbY05HdeusGvBuwiVw
8Vw+Ls/N2ti8uB5Xfe+/KN/Kh2PJ0lH1mlpyjt44q7ONjyE5gy+F+bXoeG4v
k7Wdj6MP/x63SK7Fjq5t16N38vFbUM95c7MWYlv7jC7t4UNcJSNJRa8W8SN6
3QqH+LBf7aIU/LIGNXa94p6OfPjYTHq+yK3G7ylm+TcnPjznHJIfT6rG6iiG
2xlnPqKLhNrN71fjbtXd9v2X+fjqGVS7xq0aVsu2pqp68MGMVum+rV2N7+13
trF9+VDo7lcdKa7Cov2bL6yL4ePAoxjLuNZKmM8v7ayM5SNH7fl3v7pKhFO7
zI494sMmpop5vbASohYHNP0T+OjecPpOSlwlRk/Z/X6XQuqLhUdU21SC5xLp
H5BL8nHr99v3vQJRkfSjjgby/fbSmjjpClAWV8Vdm/g4+6TVoWB2BcTlJ/+Z
2cJH4+FVIp8nyhEYP+uMJsXHhaZiTe+ecnikLVEPbiP9+bFPxOtFOY4Wmr7U
+sRHpnB0S8eBcsh0Pm26PcVHqF3nW+v/ymBrpac8Y5oA0y5eqPG8X4bs1maP
f6YLcECn7GtSSBkMqR9L3EQEYAr2VK/wKINb7TbXk2IClHHF5oyYlYH5giW9
XpqUv1aaGStaBgWNs07PZAS466IRd3m8FHaZo6UqcgLkUMNfm76UYip1sYOc
Aokf2THWX1OKNY/OMYSUBVifnur1/FYpggP/nOKuFaBNuqnbblopOMKh+WYa
AixVOP4hb6gES3yVZtevFyAxOWfuMn4JXl0zyC3REuD35jXzTlWU4OPlMJH0
zaR9GbmjjoEl2Hl82VPv3STeqcGtUqoEoW15kyOGAgw1xO5TnVmCtiOG5peM
SH+N8l6m/ynG5QNO42dMSD6tFoekPhQjfvdLU9ODAqRfHTz1T0YxxjT2Di85
ScrlDpc/NS6GXyhlWGMtgE1P4AC1rRhzvxyPu2AjgJFg6Qm1DcVQSrq4q/CM
APb5z9btky2GvtSDKHNH8r1dFg4P+orQcFlx4LeTAHWixx7O7yzC4aZUJDiT
8ZOpN3zNKsL5wMK+gcsCzNofY+jPKMLt311bAq6S/j5LvKN4rwgLj9jfWe1F
8ulXvfM9sAiPXtA9LG/S3jtFo36fImQ7ioQp3RDgZu7N6bYORaA61D+89heg
O2a56aFdRTi+OU/LOlAAz8YnbO+tRfgYqRcsEkz6I3zPqEmzCKP7zTYc/FcA
3dnbZX4sK8Licje//ntk/Hr1f+sLFeGJ4lTr3fsCFHjXJ1WMMbDOO2iNTqQA
l3WeODv/YGCHVizXL5rEH7Pw2fWJgbMppSsUE0h+U4bL/9Qz8G36Xp/KRFIu
vtODqmLAw5piOiSR7x88KdtawkCwzCevVymkvZOeNrZ5DGSGiDWaZQogG8KP
9Y1nYFPv/aWjWQJYGn6b/iqGgRJ9RfdHz0l9t+aeRZEMMMc1lL7kkvU0fXKN
ZSgDxywLXe/kk/ihxPgtwQwI8vVrN70i+djVJhj5MzDifOTyrUIyP29Lzn72
YeB6fVfVSgZpP9lX09OLAdEV9nItRWR98Z8lb/Rg4I4f7exeIsDK0YXFClcY
kO32rlhcRub/RYjf+ssMPNYTkakoJ/E7GscvOTOwKibMyb6SrKdVd1a9d2Qg
76d0mXi1APPNvKRc7RnQM0+UellD8uerlGw8z8D/AEiO9uo=
       "]]}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->True,
   PlotRange->Automatic,
   PlotRangeClipping->True], TraditionalForm]], "Output",
 CellChangeTimes->{3.52985541603708*^9}],

Cell[BoxData[
 FormBox[
  GraphicsBox[{
    {Hue[0.67, 0.6, 0.6], PointBox[CompressedData["
1:eJwVkns41QcDx106ncmQ3OulMj0Z3U4jr6KvypQ2l1xmikYajVWava8pk0WR
zqnIZS69nZN4M5fjd3PJJcfcTxyHOO4c0lihhZlIa398n+/zeb7/fZ7v5lPn
3b5WUlBQOP0+//RXdeGN0VvPoFG7zFsjuQiHO9NT+IZXsDS4tNFSWoRwxnPR
j3MTyxk2Z1KMhOgQxxfdNk5BHlwQ7yXEi8rAkdBtmfCRnRXxBEIoehlfN83k
45LZhPsa9WLwXDTv6/OzcUcyuUvzUDHWH1Go+CAnF7nbDWWOl4uRYzfT+Vde
HjKpD7Ndq4rBsR56OVFUgAxBn4aXAYFqTiurlxLCfA/rJOFKwETn5K3TPgTW
jm/QHYkhEJzYlTkwRMLBj8VvbyFAqDk99PCn0eT6g6/DPIGF63X0kzEG/f2f
2q4YkLBZbSOyDyxFPq9K/KcFiStXqNbKiTJoKa4aFviQaF4x67MIeQT6ZoWo
7TIJ9Uv3nxdMVcDiy+oTF1NIeCwYzJqEVkGvNGifaiOJjLDEd1mvq/GxuvtR
+0ISr8JXXUvYW4POirCA6DESWnvsuA1LNaiz79tjcoHCnrlLSUpVInj5fHVf
vYeCN1n68/6oWrB//uM34wEKkaGz/7uIX9HMV0vXmqdwb8fOnBKFOpyziWJ+
Z9OofRmc/1pUB+GUSPjRFhrjv+QSO2LqYXaSnFOypcH+ZrQ02L4BoSVj+etd
aJhtNarOZTUihLvTcd0ZGs7j3nWjDY3gpqUpRsTTuJCd0mIU34T9GUaf5yXT
SPaXth93bAbLyDv6VSGN0o1qstQ1LUi8k6+mX06jd/DIYIe4BSbKKnJNKY3l
zNgxdZ4Ynp/bbaqfo7HxeM3kUecn2HR6kN+xmsFB/eWZaxqt0DsvCl9azyCw
2+rP2vZWlEwZs1QtGVxPDlteSWyD2lGNxzw3BgVuQqV97hKs420beuvHICj+
P7H71rTD1zpORZvLICLoTlxvcDvOi3knigQMEhyIhHBxO6wx6+vJMMjaIuHp
bJPiptWxteYSBoWrpm5TXCm+ccwZT3/Pj8fWJB+bkkIlYXeaspRBe61p2oxT
B0aOnOqOaWUgFzhk8Io6oLOPrd/WwmA2+vRdc41ORCqVtx1qZMDyu8JvPt+J
B7sVFaZrGeiCnx3U3ok5U1EJt5yB6T+6OU8RkMJJdhYysH7bn5ed+BSLdi9O
7L/L4LOBxYIDr59i909bd9tFMfCp0CseduuCcOfg27euDM5mWFI/Ul1Q7PQL
CP0Xg6gI95IN2t2w7tp+zUefwe0vL5SXf98Nz9ZnM4vKDO5b3ar06urGgDTI
+eIoDUq38PG8pQyCckmnpIZG/XxL7Z1UGbQtffnmSTS6n07UcxZkMP5lV1p+
CI0JanWzxKsHbjruKkGHaLxJMnlytqwHn0z0n5EY0VD97qBE1aAXusns8rt/
UTA85teRF9GLuJfVDoXv/7pjV1TX4b5eJDefo2fFFOw0snrG9/ahUsXueeQD
Cq7T5f0xmX2wMD/47mgshVOtsqHNy31I0uda6fhTCCuYlz/26ce312uWb9hR
uHpDa9y3qh9Xx3tSvXQopAVzJpYMBxD2XOffh2UkHjq6vEiPGgDL2j79EUXi
kenZaavhAXx4QLZIxJEQs2/80YVBhLzhrBV9TWLw+cO5MP4gnlU0LY3Yk5iu
b1jQVByCi52nVG5OQiHn2Ruh/xA8nDbaJ7FIrItVWnGqHUJjxkpW6RyBjwI2
Kb40HsZcXJKt928ELA7uX5UQMwyB8pzAqJ/Ap5t92KbPhjHsPKxm/CsBcwVb
VoDxCNikctoXVQToGfFq0m0EB6ZqF8+VEbAdPv6BYswItkub7mWRBBraJlVc
qRFsCf9v5RcPCbhU/6B6b2wEA5HrbeoIArJCttq0lhzzVr7/nywl4Hc3Vd3W
Xg7Oj7eu/vR+n+RuWcv9Xg4/QTTlmUfgu0has/+BHO/ETs6PpATehBzSMuuS
I8rQY2G5m0DMiQ7tCNYopEXmRt/KCKh+5q/bZDkKdixdWLNIIHnvKz29wFF4
PKjQ+H0dCUOzywaBqaPgxBqVKW4jkWugvoFpGAWL6zZec4TE33WEHq8=
      "]]}, {{}, {}, 
     {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwV13k4VVsUAPDSYKgkISEVifTwUpGiFZUxFCpToolQhkSmeuYSijJFhCLj
dXHNyRQy3QGR8Z5LQshQhNJb/rrf7zv323uvvfZZa5/dVxwMrrOtWrXq9epV
q1Z+PbbnpO0o0jp+XHAizl6PBHfbO7a8j5IE3/3vpL30SbA94EXK4yhFEOlR
tAw9SwK7Pq9xpygNGCw4/K7AgATPOSd5rkddhJHXXge2XSTB28S9KRZR1jBW
fvWZtCUJ8sbU7M2i3EBO6ExW0B0SWA5NyHAb3oNJmaMCTS4k4O6P+V652R3i
s+4I8LiSwJY+6bz3kQdc+Z7amXCPBLuLX7hPe3rDzyxuPao3CRrU3Gtqcn3g
ilBXuOdDEijy0GkmUUFQZHs3/k8CCZz1nbqSfwfBGWNxhYhXJMgK28L8duUh
qJcHxEglk0Bsk8F3b7lHkKnjXWD6BufnatuU3BAMyrKFh9ozSTC0pkN7bDEU
6snjucIlOJ+n/ZV3chFwzlHZeL4d98fBsVf3eQzcEFHa1MORC2zTnVXcteiC
vluNnLnw1BnSqLMx4Gs++6CUKxey7nI7nzOMhT/PhIriN+bCkGcW+3neF2Bk
Nml9kycXjB6OyJs/jQPpAxNB4oK5IJ9k+cg2JAGGqkrMBaVywe2mVNxJtWRw
PPGYOqWVCwb0U2E3RN5ChdCXbJlnuWB3edOk390c2BXrYhD4LRca+a0qBQ7k
wZ+3ulb5p8jwJ2xL5cgNCgg/qVNuyCDDCfEu1rvQYoi5MzUfzpsHwZ/4dprc
L4NR44413S55MBufFu6UUwH2HN+3yQ/kgbxqF7dUeyXUD256mHwiH5yrrIdj
tKqhd2k2ZDk+H1R1xEwa+mtg8XbXlv2/8+H5UHRXt/EHKG7zUv2lVwBuQUdm
M2fqQEJJrW74VQEQitcSl/0bID9roFZ/pACYDA3XnuAGMKl98FdjFH1rv27x
0wb4K3srF8bQKTOLTvENsMAe5CIzXgD9m30ufMlvAKnqXT5sUwXQ+/XlpmZW
A/yIqgiKmC+AruhOz9gTH4G798oBrfUUoM5rGx/6/RFoT2pJHGIU+HE8bt9T
tkbIaZA+8R0tFPht8RtHI/C2RYp3iFPgBl/IyxT+Rjgr94mVIEGB5X9bWLxy
jcDB/+S49D4KyNrq35q2aoTbL3k1d/1LgbBeQ9+c+kYYCb4b81CFAgXirw25
Whvh2ultx02OU6Db9seeG+2N0PSeXX4fUGDvwvP6HUQjbJs7s7b+BAUq+D9t
DF1shHh6ud7iSQpM6hlH28k0QQz7EwvQpoBulXmW1LMm4Hjg9VTtAgXyfgxw
CsY2AXfRQcsptIDUVWv2xCYoYDSEv7xIgf6wm2LDGU0wwXFYaM6YArfNXWNS
qpsga3fJu0QzXO/8E3/RmSa4ayWk0mFJgVaZGnO+c80wFGP5rMGWAvJW6qVr
LjaDzdsgtyt2FIh6/nHbrHkzWK3bXbaEtliiMug2zRBnb2P5zy1cb0Ov5pP/
mmEtyW93oAMFuK/NHeLKbQZ+A3VHYRcK6L/Yt2nV5haotgjvafbC/XDffv02
fwso6FoXXvKmwDUTzne9wi1wryeDMYG+JzhqXyLVAhe6VKkbHlAgOepts7Na
C9yQluk87oP5i5AM+XK3BbJG0gR8AigQ81hiQ1NvC+gpv1znGkoBMTv+q0qD
LSDwoWduEp2lva4sbbQFvE/LSVmHUaCK84ut/1wLcJxb0D7/hAJjQa8bVXha
4YDi1adS4RRQ9hcPzj3ZCt3LMfyZzynA9NrNGZXRCtI6+RM74yjw6uFF6+Hc
VhAQOivwAG31PPSDQlEr2BX6BPShWZkLPp01reAt9loxJp4Cg930RcG+Vpj/
Oje2NoECw0d8x+M2U2HIlOtJ2SsKpJ4q1h7np0L4Hvfdm5MoYH128q2yCBX2
0KMXrdAj1qbX+6SoIPZx5r91yRQYjZLvF1Wjgva1kP6TKRQY/8GiJrlQwVzl
KvHyDca7arvsjAcVQs8MML+gb23UD1HzocJ/1aaSMqmYL/EyzcFQKhwhbB+V
oL8bPKvak0YFpUPFyR/TKDCTezI/7TMVYrfM/lufjuet3H3LrwEqfLBbTFmT
QQHnBpKD5jAVjob5WwB6dkD4n9EZKvRu1MooQHuwexjwr6dBoSblZVQmBbwU
DJou7KQB20z7H8VsCnjmBVZqitHg5fnTmdfRHrJllKMSNLhK1BRHoN337nkl
up8GZ7YK/BhDuwrMuQwr0ICRqTrxLIcCdyOkbbuUaBC9+Xx8OdqF+/LlRmUa
TFRSM4fQd9Y3aOWo0eDxHelxeRIFHOdiRV31aHD32MfZerSDc+tWm3M0+Clu
m/YNfXuSjdPUiAZJazrKuHNxv4btZlVMabDcrMxuiLb9pNKw7gYNxhs+7mtH
3zRwfjdvQwOOFNkvM2ib1tS8UTsapIZ/XL2FjPWjbvPLFica1AR0W+igrxYS
TpHeNHg9sZxejL4iL2Ad9B8NRh6mOzJWnKNt7u5Hg/LivzHf0Jap+eqXHtEg
w2ycWySPApfFRo7ph9AgdjLJ/CDaIkHkgOoTGuSfV92ojTaPChCWiKRB9T3v
8LtoM95Snm0xNIjxPmodjDYNm1zHGUcD9ZHdiQlok8CLU+OvaMAtmHf4A/qC
674PBVk0KOAl/eXKp8D5mUulqSQaPHeJHBFCG92OIMXk0aDpjOwhabThjaVY
r2IanBSvmFVHG7Dkntwuw/j7u6yNVmxxzd+yggZPgieOW6HPXmi5faqWBsGV
ugLuaH3G6msK9TTY+0+FgD9aT0/BRKoR863j7BG24kZbPaEWGrSNeqnGoHXV
E09upGG+5qfsktBnqtuOLDNooBxMLKSv+DiH7FQHDbw+W30jo3VKlcVZXbhe
tYATJWhtBSfB9h4auGhcXniPNufoDV43QIOhmZSIOrRDt/ofBRYNZlKsXJrR
vllkB5svNPCrIyfQ0c/vi7BiR2ig9yt7yyd06tkgo6ZvOD6HY8tndLHYTN3v
SRpILWxh9KIbf5gryc7QQOF39s4BdF9dfeblnzS45GpRyERPxciLhv/C/fij
EkWg2exePq1eooGEifX7FfOpcKz5sUyDRZ+fB1a8d/OduxJsdJBPXjO5Mt4R
ou/rhXV0OM6fPdu3Em++pulDDjo0W2862YO2CMhvLtlAh7yLJ3o70Y4XReEb
Nx2G26zet63Et+8RWYSXDja/Q8dbV+JbmhXX46eDohhx9eNKfC0WUQ8Ecfx6
N4kadEniRw6yMB00Em4fKEc3Ox3yZInS4VswK6gA3X8ycWKrGB18k6akslbi
GXFhuErRQTX8rE7sSjylA6fe7qeDemBf80p+JUO0iz7L0qG91yPWbyV//+6K
Vz5Mh6AlCz67lfWzPea+fYQO66+yqi+trL/953+Jx+iQavu9XB8dea/pOpsa
HTwydz8/gH6rrdB18DQdgtcFOO1Gl4okaV/XpIOu6uQLnpX1VbrKNejRwYx0
ljGO530qgkhaOEeHBZ9HzM/oNdfP8O0/Tweaq8axOrQUp9hCqBkdivjaP8Wj
j/aE2L23oMOjnsFtQegz2fN9U1Z0mPm7+MYR7XSupdrQhg5iaTXpqujS2Hsh
2+9ivnJuGPbg+9xsN7isfY8OiUbGqyrRAyp6Tl6edDAlWSyloNeyxC8M+NCh
JeNg7020rjR155swOhRspZ2ZwPrCLJXI//ctHa4ozXG0Yb3iXZPxbX8mHSaP
LBwhoU/qyO6RzMH86ndQg9GpPQqROwroQA2PkD+Btvut4cZVifv10puagvVx
TsX22FAnHZQXbXz1sL5KBU7eGeimQ+X+x7K70SatzlndfXTIOet9cjaLAuUW
njvogxhvt+HbSLTPg5C/777T4YbkFas2rO8bq3JqotkZ4LH1UI481v/jnAd/
R3AxoP/m4NwS9g+Hc0WHwjYx4PMzz4xaNIN4/8Z/KwNkZaYNjNDRbPQgp50M
yFCuzbJ/i/371KyOjiID6rg3pvhgfzpSr9i+fIMBXp6DMSrY//wlBjZ/u8mA
9vntP39if6T7Bep02jPAt9SuLhtte6KjmuTMgISNsw3C6PgSZ7LlfQaocEee
m06kwN/M7LDqSAZcsEoN931Jgbqn4lqBtQzYNdX+6XgM7u/3Rn/negZsUlPf
ORyN9VfXudKikQERf2KIEPQvzipFRRqut863vSuKAtJ+FntHehiQK9fqZhtJ
gdC7sWu0ZxnQOkzvdInA+mm6uWKTWBswdrnvlAzB+yH3hKuSRBt0vt2zmvKY
AltrGuWuS7VBqY2wjhqasT8wqVy2DS7P1QuaBWP9/PPb3+5YGzDb5ZMCHuL5
eDWm89GoDf6aOkRX+FNAffjDZ/+gNji8Mz69Cu9XR5y9fi6Nt8EsVBtux/ta
Yfwdl+nvbeB/L2xnKN7vDtfbzg7PtGH/+Mdl1cp9T8Rkmv6rDSKDM4u/2FBA
pu7wRNradojOuVKQgd8NYkLfhwxF2sFYyei16BUKbKi2as/UaYfVgrROEt4v
+7do5JtltINIxIPLfqewHimkVbxw64Alq4Ezalux/xPNEUKqn+Cn9+kAr94C
cOA7O7HhzycQkDctLHxdANLNNpG9DZ3w+0/1rO61AkjbMFAjGtQFItOfxeMk
CmCELlSWpPoZAg9rnaobzocpkfG/azZ2Q/HpzYsKyflQ+l2oqbOhG/Ltn8V9
v5IPZzQGBoMje8BMOURXYVc+vNq72o9k1AuTg9/2pgzlQXHIyBgfTx8c3K04
n5idBzFu11+Z9PTBHeFf/0S75oHSGPRuju0HK3aV/FOaeZBMkK6/uTEAHiqL
F46I5EFc0a+K5psDoDymR8QI58HzUFXBH/YDMLHA0bkglAcBSm2NancGIChZ
fN277XlgEz4nx3wwANL8Ug6a2/JAVlVlUSh2AHx9k5/dwe+vkleNYU+bBqAg
tY0pzJkHI6SZHhVBJkxqGfxTOk+GHx0pdoXbmXCkw04J0Kt+Gy3JCjPhqrF7
eu0cGQQ1i4R2iTKBtyk/j/aTDBpMTxO2PUyYFXO3m5glQyrPus46OSaY+/r5
K02R4ZrTNoa+BhMMotdUiYySwTG6wapBkwkfx410ckbI4PXOffqENhMk9vkd
PYF+ztnLI6/LBB77LI3rX8lQk/RKn8+QCZvoqaIFX8ggxtjX0mXBhIS0h43X
WGRgyis3WLky4U3ot9u8vWSY83sosOjGhIVxDpe3PWTY0NF+LcKdCV4U7+nj
aAW3W6trvJhgr550/lY3GULKEo7u8WOCZY58LbWLDIqnVmcPP2HCO6dtpzM7
yBB2vi7cLp0J98pFxN9RyfA6lZe5JpMJndr/aV1Bl8xbyMZnMcH0b8QoO3oo
dr6xmYTrkf0aZdRKBqV+qbWyhUy8X+xPm23G5zceu05VM0HuVfFrzUYyHLt3
1uxuLxOet4+f0PtABkv/9JCZPiZ0HfyYyY72f8pW4TDABFcTS5+qWjK0vi3Y
ZctiQrNWk70C+spnwWGLESYcmB88JVFDhuCjhKPmDyZIbdL8LVBFBpL60eQP
P5kw4MXJ11lJhnaDZ21q80y4EJgRFY0WtTutoLLIhKiep73b0Xlx6UsHVhEg
1cBtJ/6eDJ9/OwUKbyRgecjtmfY7MixzNBXFbCKgVzzBZAtajH/PKP9mAo67
ZAd3lZPBXuaTDg8vAYbfA2Rs0Kssjm5ZJ0iA9Si3a0gZGaTes8VPihMQv5S1
ZbSEDG4+z8hVxwgwq525MVNIhgzB+9cUVQg4Zm16tR7dS7LZln2cgIW+03Xx
6BP9Kt4xqgSImr94r4nmOjai6aBBwBCV8eANhQwJP5SZIoYEZPcvyjsXkIH2
WPJZhBEBgZ5NkbroNeK86hwXCHgp1WexD21z7mvmrDEBk1kkcVY+GftouFuj
BQENzAduxug662Hue7YE0CKn1YzzyPBrmVY1YUdAhYYm+1H0/qgyl6u3CJif
9z8ggn5a+7Rbz5GA7SejiggyGUx3H0uVcCVgus38kRM6pFjCJM6NAD1HKa8L
6Pf6PBu3uBOQ75nSfgy95/4Xp9+eBNTQgwvZ0ZPdT1TafAjgHNJVeZ1Lhl3O
HlOafgSMcnwwD0Ebcl5PqfAnwEiARLigixWPcmYEEWA8VCasgfZ9PtT+XygB
pBmG/hSJDAX7qUFzYQREpf2O7UF/rS45av+UgDv0C1r16DPTYa8uPiNAPNi9
MxHNr6dkLxtLQHfk/dzzaAOe9301LwgwEJTKUEOHMU7pm8QTwByOFvoXvf7i
WXn/RAIcJ2Ivb0Srbf+Usj2JAJnorYeWcshwv8eMn5RMQKij2IMx9Nxl61+f
3xDwWr7lXiNaXmzipkMaAZsf60uWoW8POfesTScgcubgmSz0sM39Crks3B/7
yV9P0J+cIv0D8gg47Jnfa47mPSQ8J1SA5yOOpnwOrTf3yjqXQsCZgEFedXSd
R5Z2TzEB9op5uw6g2VTkyx1Lcf8czhlLoo//LZJZX07A/XPZ60TRhX41PAfe
E2AZbVm1ET1zWsu3rpKAs7deUNeiZTmos2bVBKxlnDf8k00G20aj69M1BIST
/M7MoVNDuj8FfiAggMZX9h3N0rPUFKnH/A0sJIyiRbcMl5Ab8DzcFP81iDZt
s9uv0UhAR5dPXT86KnI6vreJAEXPdVzdaMZFN27nFjyP3gnVHWhuoT8P2KkE
vJlUnqKjtXt9p+NpBDiX0yJa0YEJHFflGQSQC7UymtDVlmHt9W0EJIZGK3xE
L4vxqV/qwPkZ6Yr16KNfYotmPuF4cw7ZH9CuaTv3PezCfOp+jq5F591882JH
N77PXB1zNejJ/fs35vcQcETNoHHF0pO53pp9BISNAe/K/2/kKnzv6ydgR3dI
64qTncst7zDx/WyQ/Lsyfv8hNQYHC5/r/01ZmV9ovv5kwiABrl9Xv1tZ34US
XcrBLwR8G9+h3YyO8Gzb+3GYgJE4VR0qulXFJMZiBPPNa1bJQHOtGuD8MUqA
v55x+ie0evU1z0ffCLA9ILm+B+3rPzYuOkGAlxfl0wC6Qt3RomCSADpjSugL
WqHJU3VgGuvb8z3LU2jn0NX5LrME7MqZj59H5+gH7eH6SUC7jAxpGb23PYL9
8C+sT/Z3ZLjxPFyJEnRvXCCgdvhSrAA6wThh7PISAR+OWdzbid7Wl94SvEyA
ZKLBE/mV/5sqi61bxYKhwJd8yuiczlbXB6tZIOTf/ec0Wp0xu9NlLT5nG1sy
RUfoB7pMrmOB/0Qczw10X7PgRxt2FlS+X37ohHapV3G+xMUC5dnnrx+hK0/R
6jo3sOAbe5pGFHpD9RVhg00s2GMnqZeCTi4PqlXnYYHtY7mACjSNTBf4V4AF
HFoPYhfQInLX7DK2sWBH4N4LnPg+W2fNvd+znQVH3+Td3Y5eThW+uV2EBab7
DvYdRf8Tf72MTYwFrLpBHV/0w8CFyx0y+PymbIw41qP2NSEF+nIs2OdV6KGE
3ukjytn4LwuiWH8r9NGFnifzKg5iPDXB9t7oIcfQtW+VWOBSY2PSh1Y1253u
cZoFf2kqN0kr9bQr/88PdRawb/1o0YTuOq9u4KDJgqIjvQVf0Y5n7Zau6rBA
itD03Y31OeE0RVf3HM7X/ScxFr0opzWz8xILuOyS2+Ow3vuFMNTrLFhgyLG1
ugq9cdQszt6SBR3W+/aOoEWTb58qucqCxaHWcwrYP9T4nkcZ2LKglCvjYBc6
+Ff/0QA3FizJJa2Xx/609bzNk/3uLJg8v+HsZXQ8eXqQ7sGCk+HS3KHoHNu1
oaL3WaBbyz4yimb07hso8meBRUoed3oRGYSrXPzGwllQLMSzQQX75esdy51P
n7GgycnzqRNa1iPoH8VI3F/S44ep6BMHX3T4xWD+X/DStpTi/evNe8kdibje
zIXiCXTWI65m/SwWJF760FaK/fnw8LNdc9ksuOWoeWcOXaG24248iQU301X9
5LG/05bkREfzWOBZv/5HJvrHrfOOviUs2No+zPe6ggzejf21UmUs4OSQPDSI
Xi9ps51azoKcrtJqMbw/CDI9qoUrWUAe0JJJQScpr91WXcWCi/xLHENo6dhQ
O5saFpgP2xlK4P0j/6dAJfcHFohyuW20RisbvOKj1LHggPD8wXT0/57/jy8=

       "]]}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->True,
   PlotRange->Automatic,
   PlotRangeClipping->True], TraditionalForm]], "Output",
 CellChangeTimes->{3.529855416038631*^9}],

Cell[BoxData[
 FormBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0.7028846900387635`", "0", "0", "0", "0"},
     {"0.0730032194381388`", "8.012451976206838`", "0.2559565886873755`", 
      RowBox[{"-", "1.4776005438794475`"}], "0.3033112282529834`"},
     {
      RowBox[{"-", "0.17386726982550355`"}], "13.664447789459537`", 
      "0.03972545781107865`", "3.3222243248532077`", "0.3331814041561628`"},
     {
      RowBox[{"-", "0.34391241976909637`"}], 
      RowBox[{"-", "10.572766082605037`"}], 
      RowBox[{"-", "0.046365175837649104`"}], 
      RowBox[{"-", "0.8606423763631201`"}], "0.6811818184219889`"},
     {
      RowBox[{"-", "0.01409271867886159`"}], "10.731927329236091`", 
      "0.002231171364957409`", "1.4823705426670182`", "0.07523792271842314`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{3.529855416045958*^9}]
}, Open  ]]
},
WindowSize->{959, 1024},
WindowMargins->{{Automatic, -9}, {Automatic, -8}},
ShowSelection->True,
FrontEndVersion->"8.0 for Linux x86 (64-bit) (November 7, 2010)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 2125, 49, 164, "Input"],
Cell[2673, 71, 6030, 146, 620, "Input"],
Cell[CellGroupData[{
Cell[8728, 221, 16413, 446, 1398, "Input"],
Cell[25144, 669, 116, 2, 31, "Output"],
Cell[25263, 673, 4198, 75, 238, "Output"],
Cell[29464, 750, 10119, 172, 224, "Output"],
Cell[39586, 924, 10975, 186, 238, "Output"],
Cell[50564, 1112, 10142, 172, 238, "Output"],
Cell[60709, 1286, 10753, 182, 229, "Output"],
Cell[71465, 1470, 1276, 29, 102, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
